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ABSTRACT 


A computer code is described which yields accurate 
solutions for a broad range of laminar, nonsimilar boundary layers, 
providing the inviscid flow field is known. The boundary layer may 
be subject to mass injection for perfect-gas, nonreacting flows. If 
no mass injection is present, the code can be used with either per- 
fect-gas or real-gas thermodynamic models. Solutions, ranging from 
two-dimensional similarity solutions to solutions for the boundary 
layer on the Space Shuttle Orbiter during reentry conditions, have 
been obtained with the code. Comparisons of these solutions, and 
others, with solutions presented in the literature; and wdth 
solutions obtained from other codes, demonstrate the accuracy of 
the present code. 
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INTRODUCTION 


Calculation of the convective heat --transfer-rate distri- 
bution for a given configuration requires information about the 
character of the vicuous boundary layer* The present work discusses 
a numerical procedure which provides solutions of a thin, Jaminar 
boundary layer bounded by a specified inviscid flow field. The 
inviscid flow field if- required as input to the present code. This 
code can be used as a tool for analysis of data from wind tunnels and 
for the extrapolation of the correlations to flight conditions. 
Furthermore, it will serve as an instructional aid for courses in 
boundary layer theory and convective heat transfer at the University 
cf Texas at Austin. 

For compressible flow problems, where the wall temperature 
is constant and the static pressure at the edge of the boundary layer 
is constant, the convective heat-transfer rate can be calculated using 
the Eckert reference temperature relations [1]. The heating rates 
calculated using this method, which does not require the comp' te 
solution of the viscc ^ equations, compare favorably with experimental 
measurements even for flows with moderate, favorable pressure 
gradients [2]. However, this technique cannot be applied for nu- 
merous problems of interest, e.g. , those involving gas injection at 
the surface, nonuniform surface temperatures, real gas effects, etc. 
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A more rigorous calculation of the heat transfer rates 
could be obtained from a solution of the viscous equations, as* 
suming the similarity conditions are satisfied [3], As noted by 
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Hayes and Probstein [4], these conditions are satisfied in the 
following cases: 

(a) for flow with a uniform free stream (i.e., constant 
pressure solutions), and a constant surface 
temperature, 

(b) near a stagnation point, if the fluid properties in 
the stagnation region are assumed approximately 
constant, and 

(c) for the limiting case of locally h>T>ersonic flow, 

Ue^ 

where 2* 

e 

For other flow problems, approximate solutions can be obtained, if 
the flow is assumed to be locally similar. Special cases involving 
mass injection can also be treated assuming similarity. One such 
example would be: 

(pv). . - 

yx 

where the injectant is the same gas as the free stream. 

The proceduies described above are still too limited to be 
of use in analyzing heat transfer data from many experimental pro- 
grams. With these requirements in mind, Bert in and Byrd [5] 
developed a finite difference scheme (NONSIMBL) to obtain solutions 
for a xaminar boundary layer for a flow which is: 
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(a) either compressible or incompressible, 

(b) possibly subject to a pressure gradient in the 
external flow, 

(c) either two-dimensional or axisymmetric , 

(d) of arbitrary free-stream gas, 

(e) possibly subjected to an arbitrary injectant, 
injection rate, and injection-rate distribution, 
and 

(f) bounded by an arbitrary wall temperature 
distribution. 

The original code was limited to a mixture of perfect-gas, nonre- 
acting flows having a constant stagnation temperature. Their 
solution technique closely follows that developed by Marvin and 
Sheaffer [6]. 

The NONSIMBL code was later modified to handle real-gas 
flows as well as perfect-gas flows for either variable* entropy or 
isentropic flow conditions at the edge of the boundary layer. It 
has been used extensively with good success but the computed 
solutions exhibited pi'obiems in certain applications, specifically 
accelerating flow past a highly cooled body [7]. The problems which 
occurred in these cases were due to the relation used to evaluate the 
transformed y-coordinate. 


OF POOR QUALITY 
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where n is the transformed y-coordlnate using the standard Lees- 
Dorodnltsyn transformation and a is a scale factor which wm 
assumed constant. This transformation allowed numerical 
integrations to be carried out over a fixed interval (zero to one) 
rather than the usual Interval in the T)-coordinate system (zero 
to infinity). Therefore, the need for an iteration to define the 
boundary layer edge was elminated. However, the coding procedure 
also eliminated the iterative procedure which required that zero 
velocity and temperature gradients exist at the edge of the bound- 
ary layer. Note that if the scale factor a is assumed a 
constant for every station, as it was in [5], there is a unique 
relation between n and n. As a result, n , is forced to 

be a constant. However, the Faulkner-Skan solution shows that, 
for a boundary layer in the presence of a given pressure gradient, 
n^dge ® constant but depends on the value of the pressure 

gradient parameter 6 C8]. For small streamwise variations of 

8 , a constant value of the scale factor la not a bad assumption, 
but for accelerating flows past highly cooled walls, such as in 
[7], it is unrealistic. 

For flow conditions such as those discussed in [7], the 
solutions generated by the NONSINBL code had relatively large 
velocity and temperature gradients at the boundary layer edge due 
to the constant a assumption. Although these gradients had only 
a small effect on the heat transfer and skin friction (which were 
of primary importance at that time) they had a significant effect 
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on the temperature and the velocity distributions. Subsequent 
applications have required that the code be capable of providing 
accurate calculations of displacement and niomentum thicknesses 
which are very sensitive to the velocity and temperature profiles. 
The present code includes an iterative procedure in vdiich a is 
varied so that the velocity gradient at the edge of the boundary 
layer is essentially zero. There is still no restriction on the 
temperature gradient at the edge of the boundary layer, but for 
most cases computed to date it is small. 

Another refinement added to the code removes the need to 
input velocity and temperatu- •? profiles at the initial station. 

The code now calculates the initial profile using a Faulkner-Skan 
solution for a perfect gas with Pr = 0.7. The initial guesses 
for f”(0) and g'(0) needed to obtain the Faulkner-Skan 
solution are calculated using a correlation of these parameters as 
functions of g(0) and B. Thus, a large range of flow conditions 
can be solved with a minimum of external calculations. 

The modified NONSIMBL code which is called NSBL, i 
coded in FORTRAN IV for the CDC 6600 at the University of Texas at 
Austin and is described herein along with a discussion of its us ;:s 


and limitations. 



NOMENCLATURE 


A, B, C, D 

*v ^ 

c 


c 


f 


c. 

1 

c . 

pi 


c 

p 

D* 


D. 

1 


F 


f 

G 

H 

H 

i 

k 

L 


matrices, equation 31 
Chapman-Rubesin factor 

local skin friction coefficient 
mass fraction for species i 
specific heat of species i 

specific heat, ^ ^i^pi 

displacement thickness including mass injection, 
equation 42 

diffusion coefficient for species i 

dimensionless streamwise component of the local 

velocity — 
e 

stream function, equation 7 
matrix, equation 35 
matrix, equation 34 
total enthalpy 

mesh point coordinate normal to the wall 

thermal conductivity 

length of the Space Shuttle Orbiter 
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H 

M. 

1 

m 

N 

n 

P 

Pr 

• 

q 

R 

r 


Sc 

St 

T 

u 

V 

X 

X 

a 

y 


Mach number 

molecular wei: it of species i 

mesh point coordinate in streamwise direction 

number of mesh points normal to the wall 

transformed y-coordinate, equation 11 

pressure 


Prandtl number. 


_E 

k 


local convective heat-transfer rate 
universal gas constant 

distance from surface of body to axis of symmetry, 
measured normal to the axis of s}rmmetry 
transformed x-coordinate, equation 6 


Schmidt number. 


pD,. 


Stanton number, ■ cVt"T ' ) 
e e p r w 


temperature 

velocity in streamwise direction 
velocity normal to the wall 
physical streamwise wetted distance from the 
stagnation point 

axial distance from the nose of the Space Shuttle 
Orbiter 

physical distance normal to the wall 
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a 

B 

6 * 

n 

^edge 

6 

e 

U 

P 

(I) 

T 

4 * 


scale factor, equation 11 

2s du 

pressure gradient paraaeter, ~ 

e 

dlsplaceaent thickness 

transforaed y-coordinate, equation 6 

that point in the boundary layer where u = 0.99u^ 

T 

non-diaensional teaperature, = — 

*te 

aoeentua thickness, equation 48 

viscosity 

density 

coluan aatrix representing the unknowns, equation 31 

local skin friction 

streaa function, equation 7 


Subscripts 

a,b,c,d,e,f aesh points, sketch 1 
e edge value 

i species i 

inj properties of the injected species 

N evaluated at boundary layer edge when used 

with matrices only 

r denotes recovery temperature 

t stagnation value 


w 


wall value 



properties of stream species 

evaluated at the wall (which is the first node) 
when used with matrices only 
injectant species 

Superscripts 

body geometry factor, k=0, for two-dimensional 
flow; and k=l for axisyimnetric flow 
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THEORETICAL ANALYSIS 


Governing Equations 

Solutions are sought for the laminar boundary layer for an 
axisynunetric or a two-dimensional configuration with possible 
ciblation or transpiration cooling. To model this flow, the present 
analysis is not restricted to the requirements for a similar, 
laminar boundary layer with mass injection. The injected gas may 
differ from the stream species and a general injection-rate distribu- 
tion may be specified. The body may be either axis 3 nnnetric or two- 
dimensional providing the radius of curvature is large in 
comparison to the boundary layer thickness, \.e., centrifugal forces 
are neglected. Approximate solutions for a three-dimensional 
boundary layer with small cross flow can be obtained using the 
axisymmetric analog [9] in which an effective radius of curvature is 
used to describe the streamline divergence. For flow with no mass 
injection at the wall, the thermodynamic properties of the free- 
stream gas may be modeled with the ideal gas relations [10] or 
with real gas properties using the thermodynamic subroutine, 

"HOLIER.” For flows with mass injection, the thermodynamic 
properties of the mixture of injectant and stream gases are ap- 
proximated with the ideal gas relations. Chemical reactions between 
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the species are not considered. The governing equations applicable 
to the flow model are as follows : 


continuity: 


3(pur^) 9(pvr^) 

3x 3y 


( 1 ) 


where k=l for axisymmetric flows and 
flow. 


k=0 for two-dimensional 


species : 


3C. 3C. 

pu + pv ^ 


3y 3y 


^ P D, 


3C, 
2 

3y 


( 2 ) 


x-component of momentum: 


3u 


3u 





(3) 


y-component of momentum: 

= 0 (4) 

which represents the standard boundary- layer assumption regarding 
the pressure gradients noimial to the wall. 


3y 


energy: 



The governing equations which describe the nonsimilar, possibly 
compressible, flow in physical coordinates are nonlinear, partial 
differential equations. Therefore, a transformation is sought to 
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simplify the solution procedures. Using the standard Lee- 
Dorodnitsyn coordinate transformation CH]: 




2k. 

e e e 


( 6 ) 


W p p . 
n = -^ — 'iy 

' o e 




where the stream function is defined as : 


= (2s)°*^ f 


(7) 


one obtains the following equations: 


species: 


[h <) » f < = (c= f' - c;j 


( 8 ) 


X' jmentum: 


ff 

energy: 
C 


+ (Cf ) + - (f ) J = 2s I f (f )® - f®f J(9) 


* ^le C f \ 1 I 

f0+ — -^(C -C 0C, + - 
r So \ p^/ 1 j 


C T. 
p te 


P 


Bf 

e e 
pc 

p te 


^ S * S 

= 2s (f 0 - e f®) 


( 10 ) 
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where superscript prime denotes differentiation with respect to ti 
and superscript s, with respect to s. Having introduced the 
stream function, which automatically satisfies the overall 
continuity equation, the new set of governing equations contains 
one less dependent variable. Except for a limited number of 
injection-rate distributions, it is not possible to reduce the 
governing equations to a form analog :>us for similar solutions, as 
discussed in the Introduction. An additional coordinate transforma- 
tion can be made, as suggested in [123 and as mentioned in the 
Introduction: 


n 


1 - e 


-an 


( 11 ) 


This transformation is for numerical purposes. Numerical integra- 
tions can now be carried out over a fixed interval (zero to one) 
rather than the usual interval in the n-coordinate system (zero to 
infinity). This coordinate system eliminates the need for an 
iteration to define the boundary layer edge. Note that in the 
present approach it is assumed that the edge of the viscous boundary 
layer, the edge of the thermal boundary layer, and the edge of the 
species concentration layer, all occur at the same n. 

Also, the transformation affects nodal point spacing in 
the physical-coordinate plane. Points which are evenly spaced with 
respect to the n-coordinate are not evenly spaced in physical space. 
Spacing of the y-coordinates of the nodal locations varies with 
position, such that Ay increases with distance from the wall. 
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This results In more nodal points in the region near the wall, where 

gradients are large, and fewer points in the region away from the 

wall, where gradients are smaller. The scale factor a will be 

treated as a constcmt for any two adjacent streamwise stations. 

The resultant governing equations in the transformed system with 
u * 

F = — replacing f , are as follows, 
e 

species : 


(fe)„ ‘"i" * ■ 'i”) 


+ a(l-n) f C 


In 




= 2s |C- F - O(l-n) F C, 
'Is s In 


) 


( 12 ) 


momentum: 


fF a(l-n) + a*(l-n)* C F + Co*(l-n) f (1-n) F - F ) 
n n n l nn n y 

~ - F* = 2s fFF - f F a(l-n)) 
p y s s n J 


+ 6 


energy: 


(13) 


— ot"(l-n)" (^) + ^|j:ja*(l-n) ((l-n)0„„ - 0^j + a(l-n) f 0^ 


C / 

— It (c - c )a*(l-n)2 0 C 
C \^1 P 2 / " 


Cu 2 

In + f*n a^(l-n)" 

In ^ n 

p te 


n P 

0 ~ F = 2s 1 F0^ - o(l 

C T, P Is 

p te 


(«s - Vs) 


(14) 
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where the subscripts n and s denote differentiation with respect 
to r. and s , respectively. 


Boundary Conditions 

In the previous secticMi, the governing equations were 
written in terms of three dimensionless, dependent variables, F, 
C^, and 0. Since the surface temperature and the inviscid flow 
field are known a priori, values for F and 0 are immediately 
determined at both boundaries. 

At the wall, i.e., n = 0 


F = 0, 0 = T /T 

’ w te 


(15) 


At the boundary layer edge, i.e., n = 1 


F = 1, 0 = T /r 
e te 


(16) 


Also, it is possible to readily define the species concentrations 
at n = 1 

= 1, Cj = 0 (17) 


The values of the concentration fractions at the wall, however, 
are not explicitly known. The appropriate wall boundary condition 
is obtained from the conservation of species equation: 


<“''>« '=!» = ”l2 



(18) 
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Physically, this boundary condition specifies the balance between 
the convective flux away from the wall and the diffusive flux 
toward the surface for the free-stream species. Written for the 
n,s coordinate system in terms of nondimens ional parameters, this 
boundary condition becomes 


C 


Inw 


(pv) C- Sc ySs 
w Iw w ' 

V 

a p y u r 
w w e 


( 19 ) 


In addition to the dependent variables 0, and F, the govern- 

ing equations also contain the stream function, f. By definition: 


f = F 



F 

ad-n) 


dn + f(0) 


( 20 ) 


From the continuity equation it follows that 


f(0) = - [ (pv) r^ dx (21) 

In summary, the boundary conditions for and f(0) at 

n = 0 are determined by substituting the specified injection-rate 
distribution into the appropriate equations. The other boundary 
conditions are determined from the usual physical constraints. A 
complete listing of the boundary conditions is as follows: 
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At the wall, n = 0: 


F = 0 



(pv) C- Sc T/ts 
^ 11 ^ 


w Iw w 

'In ■ k 

a p y u r 
w w e 


f(0) = 


= _ _i_ f (pv) 

V2T Jo 


k , 
r dx 


At the boundary layer edge, n = 1: 


F = 1 



C 


1 


= 1 


C 


2 


= 0 


Formulation of the Solution Algorithm 
The governing partial differential equations, and the 
appropriate boundary conditions have been developed. However, to 
obtain a solution it will be necessary to use the finite difference 
form of the equations. These finite difference approximations are 
computed using the nodal scheme shown in Sketch 1. 
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Sketch 1. - 
Nodal Pattern 


i+1 

i 

i-1 



m m+1 

The finite difference approximations for the partial derivatives 
of the dependent variables are functions of their values at the 
streamwise coordinate, m+1, where the solutions are to be 
calculated; and at the station, m, where the flow field is 
already known. 


G represents a particular dependent variable 
G = G 


= Si" ^®e ■ 


( 22 ) 


G = - 4 - (G - G + G, - G.) 

n 4An c a f d 


G = — i (G - 2G, t G t G, - 2G_ t G.) 

"" 2«n)’ G "b a f e d 


The approximations for the themiodynamic properties are evalviated 
only in terms of their values at streamwise coordinate, m: 


R represents a thermodynamic property 


R = ^ (R - R,) 
n 2An c a 


( 23 ) 


R = 
nn 


(An) 


- 

2 C 


- * 'a* 




M 


25 


Difference approximations for partial derivatives involving products 
are defined so that the dependent variable, at the m + 1 stream- 
wise coordinate appear linearly in the resulting equation. 


1 


G M = 
n n 


8(An) 


“ C(G^ • G >(” 4 = • + (G. - G,)(M - M )] 

2 card rdca 


S - S 


(24) 


= "n = S 4E <"c - "a * "f - "d> 
1 


.2 « 


4(An) 


— (G - G )(G^ “ Gj) 
2 c a r d 


Substitution of the appropriate approximations into the momentum 
equation, using 6i = a(l-n), results in the following: 


* 17^ f + F, - F.) + a f (F - F t F- - F.) 

n 4An c a f d 4An c a f d 


+ c 

a^ (F - 2F. 

+ F 

+ F, 


2(An)2 ^ 

a 


- c 

a a (F - F + 

4An c a 

"f - 

"d> 



F - 


+ 3 

r - '■b'e) = '■b 

e 

As 

b 


(25) 


a (2s) f 77“ (F - F + F. - F.) 
s 4An c a f d 


OF FOOR ttOAlJW 



26 


Tn the above equation, f is evaluated using the relations for the 

s 

thermodynamic properties. Although the derivatives of the properties 
are not represented in finite difference forms, it is understood 
that the values are to be obtained in this manner. The finite 
difference form of the momentum equation takes on the form 


a f . <^n . C 


^f 425T ^ ~42Sr 


2(An)' 


got C 

4An 4An 


+ F 


2a^ C 
2(An) 


2s F 

T-^h-sr- 


.r , 5 f , a’ c , aa c “ ‘’s 

" 4An 4An v 2 4An 4An 


(26) 


5 f (F - F ) 

c a 

4Sn 


g^ C (F - 
n c 

4An 



g^ C(F - 2F, + F ) 

C D d 

2(An)' 


gg c(F 
+ 42^ 





2s Fj 
As 


2s 5 f (F - F ) 
sc a 

uAn 


The quantities enclosed in brackets depend only upon properties 
and dependent variables evaluated at the streamwise coordinate m, 
where the solution is already known. The equation can be rewritten: 
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Conservation of momentum: 


^11 ^ ^e ^ *^11 ■ °1 


Conservation of energy: 


^21 ^ ®21 ^e ^ ^21 ■*■ ^22 ®d * ®22 ®e ^ ^22 ®f * ^23 ^Id 

*^23 ^If " ^2 ( 28 ) 


Conservation of species: 


^33 ^Id ®33 *^le ^33 ^If ~ ^3 


The following matrix equation can be used to represent these three 
equations, where the A, B, and C matrix elements are defined in 
Appendix A* 




21 "22 23 

0 0 






®d 

+ 


^Id 



Bii 0 0 


21 22 


° ° ®33 ^le 


11 

0 

0 

21 

^22 

^23 


0 

"33 
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28 


Replace 





■^1 

0 

0 

111 

®d 

A = 


^2 

^23 


^Id 


0 

0 

^33 

4 


Noting that d cox*i?esponds to i-1, e to i, and f to i+1, 
the matrix equation can be rewritten in more compact notation: 


A, 


-x,m -i,m+l -i,m -vitljm+l 


D. 

^i,m 


(31) 


The subscripts i,m refer to the nodal point at which the matrix 
is evaluated. The tertns in the A, B, C and D matrices depend 
only upon quantities at station m. Since these quantities are 
known^ the A, B, C and D matrices can be evaluated. As the value 
i varies from 1 at the wall to N at the boundary- layer edge, 

N-2 matrix equations of the form shown above are produced. The 

equation is not valid at i=l or i=N because the values of the 

elements in the A, B, C and D matrices for any point n,m also 

depend upon pix>perties at points i+l,m and i-l,m. At the wall 

and at the boundary layer edge, the boundary condi*^ions will be used 
to evaluate the flow variables. In general, the set of matrix 
equations for all N points can be represented using still another 
matrix representation: 

ORIGIN ALPA^ 
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Si 


0 

0 

.0 


0 0 


0 


A. 


B. 

-.1 


0 


“1 ' 



0 


~2 


?2 

Si - ° 




D. 

6(-i Sn-i Sn-1 


i^N-1 


Sn-1 

■ ^ !n - 


•!?n - 


•?N - 


(32) 


To solve this set of lineeu? matrix equations, the Gauss- Jordan 
method of elimination is used as described in [13], Briefly, 


where 


and 


0). = -H. w. . + G, 1 < i < N-1 

~1 -1 ~l+x -1 — — 


(33) 


H. = (B. 
-1 -1 


6i «i-i>~‘Si 


(34) 



A. H. ,)-‘ 
~i-l 


(D. 

-1 


-i Si-1^ 2 < i < N-1 (35) 


Equations (34) and (35) do not provide the values of 
and H, which are necessary to compute w, using equation (33). 
These values are obtained using the boundary conditions at the wall. 
Recall that, at the wall; 


F = 0 
0 = 9 
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'In 


V2s C, (pv) Sc 

Iw w w 

j_ 

a p y u r 
w w e 


Following the method outlined in [6], one can write 

C.(i=2) = C,(i=l) + An[C.(i=l)l + (An)HC-(i-Dl 
1 1 in 1 nn 

where Cj^(i=2) is the value of at the i=2 node in the 

direction normal to the surface. The value of [C,(i+D] can be 

1 n 

determined using equation (19). The value of CC,(i=l'] can be 

1 TUI 

obtained, after some manipulation, from the species equation for 
i=l. The above equation results in the relation: 


Cj^(i=l) 


C^(i=2) 
1 + ENS 


(36) 


where 





yis Sc (pv) 

w w 


a(pp) u r 
w e 


Using this expression for the boundary condition, the appropriate 
values for and are: 
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0 


»1 


0 


0 


I 



0 0 

0 0 

1 

^ ■ 1 + ENS 


(37) 


(38) 


The solution for the o>-vector at each of the N-stations 
of the boundary layer for a given streamwise coordinate, m+1 , 
begins at the outer edge of the boundary layer. Using the known 
boundary conditions: 


F = 1, e = 6 , and C = 1 
* e 1 


one finds that: 


b) 


1 ] 
0 


(39) 


The solution proceeds inward, i.e., toward the wall, applying 
equation (33) sequentially to determine t'..e velocity, the tempera- 
ture, and the species profiles. Having obtained numerical values 
for these parameters, a convergence check is made. The convergence 
criteria is satisfied if consecutive solutions of the boundary 
layer have the slope of the velocity function at the wall within 
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0.05% of each other. For reference, this tern is titled "CONVG” 
in the Fortran listing of the program. Should the convergence 
criteria not be satisfied, the elements of the A, B, C and D 
matrices are reevaluated using the average of the newly calculated 
value of F (or 9 or and its value at the previous streamwise 

station. Tne iteration process ntinues until consecutive values 
of the velocity gradient agree. Once the convergence criteria is 
satisfied, the y-gradient of velocity at the edge of the boundary 
layer is checked to see if it is within a specified tolerance. Only 
if the value of the edge shear, which is defined as 





is greater than zero and less than 0.2, is the solution considered 
to be the desired one. Otherwise, a is changed by 5% and a linear 
extrapolation method is used to recalculate the boundary layer until 
the edge shear criteria is met. In addition to the viscous profiles, 
numerical values are obtained for the heat-transfer rate, the skin- 
friction, the displacement thickness, etc. The procedure is then 
repeated to obtain the boundary layer profiles at the next streamwise 
station, until solutions are obtained at all stations for which the 
inviscid flow has been specified. 

The solution of the initial station is obtained using a 
Faulkner-Skan formulation for a compressible, laminary bound^u?y layer 
for a perfect gas with Pr = 0.7. This initial profile is calculated 
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in the NSBL code by the subroutine entitled PZ6YBAK. A linear 
double interpolation routine is used to calculate reasonable first 
guesses for f*'(0) and g'(0) based on the values of g(0) and 
0, as determined from the boundeucy conditions. 

The numerical routine described herein computes the 
viscous-layer profile of the velocity function F, the temperature 
function 6, and the stream species concentration C^. Values are 
also obtained for the following parameters: 


CF: the local skin-friction coefficient 


^ 0.5p u* 

e e 


(40) 


DELST: the displacement thickness* feet 


e e 


DSTAR: a thickness which represents the free-stream mass 

flow entrained in the boundary layer with injection 


D* = 6* + 


I 


X (pv) r 
w 


P U 
e e 


dx 


(42) 


This equation corresponds to the definition of 
Hayasi [14]. 


HCOEF: the local heat -transfer coefficient, Btu/ft* sec ®R 


T - T 
te w 


(43) 


HRATIO: the ratio of the local heat -transfer coefficient to 

the reference heat-transfer coefficient 


HRATIO = 


ref 


(44) 


QDOT: the local convective heat -transfer rate, Btu/ft* sec 


• i ai\ ^ Pw 


k p u r a(l-n) 

•' e te 


V2i 


(S). 


(45) 


STNO: the Stanton number 


St = 


p u C (T - T ) 
e e p r w 


(46) 


TAU: the local skin friction, lbs /ft ^ 

. L 3u\ <’w "e (W\ 

- r " — iTs Nw 


THMOM: the momentum thickness, feet 

ry \ 


« = I, - i-) ’‘y 


pur 
e e 


Jq “^l-n; 


dn 


TREC: the recovery temperature, ®R 


(47) 


(48) 


= r(T.^ - T^) + T 
r te e e 


(49) 
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DISCUSSION OF RESULTS 


The NSBL code has been used to generate numerical solutions 
for a variety of problems. Presented herein are the solutions for 
some representative flows. Specifically, there are two flows where 
the similarity conditions are satisfied and two flows where they are 
not. The two flows where the simileu:*ity conditions are satisfied 
are: (1) incompressible flow over a flat plate (the Blasius problem); 
euid (2) supersonic flow past a sharp cone. The two flows which do 
not satisfy the similarity conditions are reentry flight conditions 
of the Space Shuttle Orbiter at (1) a relatively low Hach number 
(M = 9.49); and (2) a very high Mach number (M = 29.86). 

Incompressible, Laminar Boundary Layer 
on a Flat Plate 

A solution of the inccrnipressible, laminar boundary layer 
was obtained with the NSBL code for flow past a flat plate at zero 
angle of attack. The free-stream velocity was 167 ft/sec; the stag- 
nation pressure was 1.0 atmosphere. Both the stagnation temperature 
and the wall temperature were 540®R. The skin-friction coefficient, 
the displacement thickness, and the momentum thickness as calculated 
using the NSBL code are compared in Tcible 1 with the results reported 
by Schlichting Cl5], and the values calculated using the BLAST code 
[3]. The values obtained using the NSBL code are in close agreement 
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with the values from the other solutions. The velocity profile 
computed using the NSBL code is presented in Figure 1. Included 
for comparison are the velocity profile presented in [153 and that 
obtained using the BLAST code. The NSBL profile is that calculated 
for X = 1.0 ft, i.<»., for a distance of 1.0 ft from the leading 
edge. It should be noted that the y- transformation used in [153 
differs by a factor of i2 from that \ised in the two codes 
developed at the University of Texas and in the earlier reference 
of Blasius [16]. However, the solution given in [15] has been con- 
verted to the present definition for n- The velocity profiles 
calculated by the three methods are in excellent agreement. 

Laminar boundary Layer for Supersonic Flow 
Past a Sharp Cone 

The NSBL code was also used to obtain a solution for the 
laminar boundary layer on a sharp cone which had a half angle of 12®. 
The cone was at zero angle of attack in a supersonic stream. The 
free-stream Mach number was such that the flow downstream of the 
shock was supersonic and the flow properties were constant along the 
surface of the cone. At the edge of the boundary layer the static 
pressure was 0.05 atmosphere, the static temperature was 1000®R, 
and the Mach number was 4.84. The wall temperature of the cone was 
540®R. The velocity profile at a distance of 1.0 ft from the apex 
of the cone as calculated with the NSBL code is compared in Figure 2 
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with that computed using the BLAST code. The temperature profiles 
obtained from the two codes are compared in Figure 3. Note that the 
maximum temperature in the boundary layer is greater than either 
the wall temperature or the static temperature at the edge of the 
boundary layer. This region of elevated temperature within the 
boundary layer is due to the effect of viscous dissipation. Once 
again, the agreement between the different calculation procedures 
is excellent. The convective heat-transfer-rate distributions as 
calculated by the two codes are presented in Figure 4. Included 
for comparison is the distribution calculated using Eckert *s reference 
temperature technique (as discussed briefly in the Introduction). 

The heat-transfer rates were nondimensionalized with the theoretical 
heat-transfer rate to the stagnation point of a 1.0 ft radius sphere 
under the same flow conditions as calculated using the theory of 
Detra, Kemp, and Riddell [17]. The agreement between these three 
solutions is very good. 

Compressible, Laminar Boundary Layers 
for Reentry Flow Conditions 

As was mentioned in the Introduction, it was desired to 
write a code that could accurately solve for the boundary layer on 
a highly cooled wall in an accelerating, high-speed flow. The 
NSBL code was used to solve for the boundary layer on the Space 
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Shuttle Orbiter for two flight conditions from the reentiy 
trajectory. These conditions included a relatively low Mach number 
(M = 9.49) and a very high Mach number (M = 29.86). The inviscid 
flow fields for these two flight conditions, which were provided by 
Dr. W. D. Goodrich of NASA Johnson Space Flight Center [18], are 
presented in Table 2. Additional points were added to those in 
[18] to improve the accuracy of the NSBL solution. Velocity and 
temperature profiler, displacement-thickness distributions, heat- 
transfer-rate distributions, and skin-friction coefficient distri- 
butions for these two flight cases, as calculated by the NSBL code, 
are presented in Figures 5 through 14. Also shown for comparison 
purposes are the solutions obtained using the BLIMP code [19], as 
provided by JSC, 

For the first flight condition, the free-stream Mach 
number was 9 . 49 and the angle of attack was 30 . 8 degrees . The 
conditions at the stagnation point and the streamwise distributions 
for the static pressure, the wall temperature, the entropy at the 
edge of the boundary layer, and the radius of the equivalent body 
of revolution are presented in Table 2(a). These conditions were 
used as the input boundary conditions for the NSBL code and for the 
BLIMP code. Solutions using the NSBL code were obtained for two 
thermodynamic flow models: (1) for air that behaves as a thermally 

perfect (p = pRT) and a calorically perfect “ constant) gas; 
and (2) for air that behaves as a real gas. For the perfect-gas 
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model the viscosity was computed using Sutherland's equation. 


P •- 


,1.5 


2.270 


T + 198.6 


X 10"® 


Ibf sec 
ft^ 


Furthermore, the Prandtl number was assumed a constant, Pr = 0.70. 

The thermal conductivity was calculated using the definition of 

P C 

the Prandtl number, k = — . For the real-gas model, thermo- 

dynamic properties were calculated using the tabulated values of 

[20]. The viscosity, thermal conductivity, and Prandtl number 

were calculated using linear, double-interpolation of the properties 

tabulated by Hansen [21]. The specific heat was calculated using 

the definition of the Prandtl number. In the BLIMP solutions the 

real-gas model for air was used both for the thermodynamic properties 

and for the transport properties. 

The velocity profiles as calculated by the NSBL code are 

compared to the BLIMP calculations in Figure 5. Figure 5(a) 

presents the velocity profiles at x = O.IOL and Figure 5(b) 

a 

presents the profiles at x = 0.498L. At both streamwise loca- 

3L 

tions, the agreement between the NSBL real-gas calculations and the 
BLIMP calculations is very ^ood. However, the solutions obtained 
for the perfect-gas model differ substantially from those obtained 
for the real-gas mode_. The temperature profiles at these same two 
streamwise locations are compared in Figure 6. Again, the agreement 
between the NSBL real-gas and the BLIMP calculations is very good. 
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The NSBL perfect-gas calculations differ significantly, especially 
at the downstream location. Note that, at x = 0.498L, T = 0.385T 

ol C 

at the wall for the real-gas solution, whereas tb<^ corresponding 

temperature if 0.685T for the perfect-gas solution. Since the wall 

temperature is specified as 1500®F for both gas models (see Table 

2a), the differenc in the two dimensionless values of T /T indi- 
’ we 

cates a difference between the values of T . The same inviscid flow 

e 

field, as defined by the pressure ratio (Pg/p^2^» local entropy 

(s /R), and the streamline divergence metric (RDS), was used for both 
e 

gas property models. Therefore, it is clear that the thermodynamic 
model used in the calculations can have a significant effect on 
the resultant solutions. 

The displacement -thickness distributions are compared in 
Figure 7. The BLIMP and the NC5L real-gas distributions are in 
good agreement. The displacement-thickness distribution calculated 
using the NSBL perfect-gas model compai’es reasonably well with the 
other two. 

The heat-transfer-rate distributions are compared in 
Figure 8. The local neat-transfer rates were divided by the theo- 
retical heat-transfer rate to the stcignation point of a 1.0 ft radius 
sphere as calculated using the theory in [17], The NSBL real-gas 
distribution compares very well with the BLIMP results. However, 
the values for the NSBL perfect -gas solution are substantially lower 
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than those for the other two models. The 8kin-fr:’.ction-coefficient 
distributions are compared in Figure 9. The agreement between the 
NSBL real-gas solution and the BLIMP solution is very good. Although 
the values from the NSBL perfect-gas solution are below those of the 
other two models, the differences are not as great as they are for 
the heat- transfer-rate distributions. In summary, the NSBL real-gas 
calculations and the BLIMP calculations compare very favorably for 
this flight case. The NSBL perfect-gas calculations, however, 
differ significantly from those of the other two models and demon- 
strate the error that could be introduced by using perfect-gas 
assumptions for a real-gas situation. 

Foi' the second flight condition, the free-stream Mach 
number was ?9.86 and the angle of attack was 41.4°. The input 
boundary conditions for the NSBL code and the BLIMP are pre- 
sented in Table 2(b). The NSBL code was used to . solution 

for the real-gas model only. The velocity profiles at s = 0.107L 
and at X = 0.*+44L as calculated by the NSBL code are compared to 
the BLIMP solutions in Figures lOva) and 10(b), respectively. The 
agreement between the profiles calculated by the two codes is 
excellent at both streamwise locations. The temperature profiles at 
the same two streamwise locations are compared in Figui'es 11(a) and 
11(b). Again, the agreement between the solutions provided by the 
two codes is excellent. It is interesting to note that both codes 
indicate an inflection point in the temperature profiles at the two 



42 


strcansKise locations. In both cases the physical temperature in 
the boundary layer in the region of th^i inflection point was on 
the order of 4500 to 5000°R. In this temperature range^ the oxygen 
molecules begin to dissociate at the static pressure of this flow 
(which is approximately 0.01 atm). In this temperature range, the 
thermal conductivity changes rapidly with temperatuz^. Further- 
more, the value of the thermal conductivity at a specific temperature 
also depends on the static pressure. 

The displacement-thickness, the heat-transfer-rate (non- 
dimensional ized as before), and the skin-friction-coefficient 
distributions as calculated by the two codes are compared in 
Figures 12, 13 and 14, respectively. Once again, the real-gas 
calculations of the NSBL code compare very well with those of the 
BLIMP code. 
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CONCLUDING REMARKS 


The agreement between the sample calculations made using 
the NSBL code and the results available in the literature, along 
with the calculations of other codes, demonstrates the validity of 
the solution algorithm of the NSBL code. It should be noted, 
however, that the equations employed in the code incorporate the 
assumptions typically applied to thin boundary layers. Furthermore, 
it is assumed that the viscous boundary layer, the thermal boundary 
layer, and the species concentration all have the same thickness 
Thus, because the iterative solution procedure requires that only the 
shear be negliglbe at the edge of the boundary layer, the code should 
be applied with caution to flows irtiere the Prandtl number and the 
Schmidt number are significantly different than unity. It should 
also noted that the NSBL code is applicable only to laminar 
boundary layers. Therefore, comparisons between experimental results 
and the solutions obtained with the code are valid only when the 
physical boxindai'y layer actually is laminar. 
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Table 1. - The laminar boundary layer for Incoopresslble 
flow over a flat plate. 




BLAST 

Schlichting 


NSBL 

(Ref. 3) 

(Ref. 15) 

^f 

0.6816 

0.6644 

0.6641 




jRe 

“ X 


1.7206 

1.7963 

1.7208 

X 

^ Re 

X 



e 

0.6541 

0.6633 

0.664 
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Table 2. > The input boundary conditions for Space Shuttle Orbiter 
night Design Trajectory. 

(a) = 9.49; altitude = 162,000 ft; angle of attack = 30.83^; 

= 215.8 Ibf/ft^; T^^ = 5508»R. 


M 

x(ft) 

'■e'Pta 

T (OR) 

V# 

S /R 
e 

RDS(ft) 

1 

0.8325 

0.9164 

2265. 

37.24 

0.8225 

2 

1.0113 

0.8900 

2251. 

37.23 

0.99C0 

3 

1.1900 

0.8595 

2238. 

37.22 

1.1500 

4 

1.5188 

0. 8120 

2211. 

37.19 

1.4600 

5 

1.8475 

0.7687 

2182. 

37.15 

1.7500 

6 

2.1450 

0.7310 

2154. 

37.10 

2.0100 

7 

2.4425 

0.6992 

2124. 

37.07 

2.2500 

8 

2.7088 

0.6710 

2097. 

37.06 

2.5000 

9 

2.9750 

0.6475 

2071. 

37.05 

2.7100 

10 

3,2875 

0.6190 

2041. 

37.04 

2.9600 

11 

3.6000 

0.5930 

2012. 

37.03 

3.1800 

12 

4.1625 

0.5540 

1964. 

36.99 

3.6100 

13 

4,7250 

0.5227 

1923. 

36.95 

4.0230 

14 

5.2875 

0.4990 

1887. 

36.91 

4.4250 

15 

5.8500 

0.4780 

1859. 

36.88 

4.8000 

16 

6.3875 

0.4600 

1834. 

36.85 

5.1800 

17 

6.9250 

0.4454 

1815. 

36.82 

5.5480 

18 

7.4750 

0.4310 

1799. 

36.79 

5.9000 

19 

8.0250 

0.4182 

1784. 

36.77 

6.2600 

20 

8.6000 

0.4060 

1770. 

36.74 

6.6300 

21 

9.1000 

0.3968 

1760. 

36.72 

6.9450 

22 

S.6500 

0.3870 

1750. 

36.69 

7.2800 

23 

10.2000 

0.3770 

1740. 

36.68 

7.6200 

24 

10.7380 

0.3690 

1733. 

36.65 

7.9450 

25 

11.275 

0.3610 

1727. 

36.63 

8.2700 
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Table 2. - (a) Continued. 


N 

x(ft) 


T (OR) 

w 

S /R 
e 

RDS(ft) 

26 

11.813 

0.3540 

1723. 

36.60 

8.5890 

27 

12.350 

0.3477 

1720. 

36.58 

8.9075 

28 

12.888 

0.3415 

1717. 

36.56 

9.2160 

29 

13.425 

0.3360 

1715. 

36.53 

9.5250 

30 

13.963 

0.3300 

1712. 

36.51 

9.8340 

31 

14.500 

0.3255 

1710. 

36.49 

10.143 

32 

15.575 

0.3170 

1706. 

36.44 

10.738 

33 

16.650 

0.3086 

1701. 

36.40 

11.333 

34 

17.713 

0.3020 

1691. 

36.36 

11.925 

35 

18.775 

0.2955 

1675. 

36.31 

12.518 

36 

20.388 

0.2880 

1612. 

36.25 

13.363 

37 

22.000 

0.2819 

1557. 

36.19 

14.208 

38 

24.680 

0.2770 

1537. 

36.07 

15.599 

39 

27.360 

0.2751 

1530. 

36.96 

16.990 

40 

30.030 

2.2760 

1526. 

35.85 

18.399 

41 

32.700 

0.2772 

1522. 

35.76 

19.808 

42 

35.385 

0. 2780 

1518. 

35.68 

21.174 

43 

38.070 

0.2795 

1515. 

35.62 

22.540 

44 

40.735 

0.2800 

1512. 

35.56 

23.920 

45 

43.400 

0.2810 

1510. 

35.50 

25.300 

46 

46.085 

0.2820 

1507. 

35.44 

26.695 

47 

48.770 

0.2835 

1505. 

35.39 

28.090 

48 

51.435 

0.2843 

1502. 

35.34 

29.480 

49 

54.100 

0.2850 

1500. 

35.29 

30.870 

50 

56.785 

0.2860 

1493. 

35.24 

32.250 



Table 2. - Continued 

~ 29.86; altitude = 246,000 ft; angle of attack = 41.4®; 
Ibf/ft^; T^2 = 10,798®R. 


M 

x(ft) 

Pe/Pt2 

T^(®R) 

sjn 

e 

RDS(ft) 

1 

1.9500 

0.8950 

3095. 

54.13 

1.8900 

2 

2.0470 

0.r890 

3076. 

54.09 

1.9800 

3 

2.1430 

0.8840 

3061. 

54.03 

2.0600 

4 

2.2400 

0.8780 

3048. 

53.97 

2.1500 

5 

2.3370 

0.8740 

3037. 

53.92 

2.2300 

6 

2.5300 

0.8640 

3016. 

53.81 

2.4100 

7 

2.7230 

0.8540 

2996. 

53.70 

2.5900 

8 

2.9200 

0.8440 

2977. 

53.60 

2.7700 

9 

3.1100 

0.8355 

2960. 

53.50 

2.9300 

10 

3.3000 

0. 8260 

2942. 

53.39 

3.1000 

11 

3.4900 

0.8170 

2927. 

53.29 

3.2700 

12 

3.6800 

0.8090 

2912. 

53.20 

3.4300 

13 

3.8700 

0.8000 

2897. 

53.10 

3.6000 

14 

4.0600 

0. 7920 

2883. 

53.01 

3.7500 

15 

4.2500 

0.7837 

2869. 

52.92 

3.9100 

16 

4.6170 

0.7680 

2839. 

52.76 

4.2300 

17 

4.9830 

0.7530 

2812. 

52.61 

4.5300 

18 

5.3500 

0.7390 

2784. 

52.48 

4.8400 

19 

5.7170 

0.7250 

2758. 

52.36 

5.1500 

20 

6.0830 

0.7125 

2732. 

52.26 

5.4300 

21 

6.4500 

0.7006 

2707. 

52.19 

5.7300 

22 

6.8100 

0.6900 

2684. 

52.11 

6.0200 

23 

7.1700 

0.6790 

2660. 

52.05 

6.3200 

24 

7.5300 

0.6677 

2637. 

51.99 

6.5900 

25 

7.8970 

0.6590 

2614. 

51.93 

6.8800 
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Table 2. - (b) Conclusion. 


M 

x(ft) 

Pe'Pt2 

T (OR) 
w 

S /R 
e 

RDS(ft) 

26 

8.2630 

0.6500 

2595. 

51.87 

7.1700 

27 

8.6300 

0.6411 

2577. 

51.82 

7.4300 

28 

9.1700 

0.6290 

2552. 

51.75 

7.8300 

29 

9.7100 

0.6175 

2535. 

51.68 

8.2500 

30 

10.255 

0.6080 

2522. 

51.61 

8.6500 

31 

10.800 

0.5990 

2511. 

51.54 

9.0300 

32 

11.350 

0.5900 

2499. 

51.46 

9.4600 

33 

11.900 

0.5815 

2488. 

51.37 

9.8400 

34 

12.400 

0.5740 

2479. 

51.30 

10.230 

35 

12.900 

0.5667 

2470. 

51.24 

10.600 

36 

13.450 

0.5603 

2460. 

51.18 

11.100 

37 

14.000 

0.5540 

2450. 

51.13 

11.600 

38 

15.050 

0.5425 

2433. 

51.05 

12.250 

39 

16.100 

0.5313 

2419. 

51.02 

12.900 

40 

17.200 

0.5210 

2403. 

50.95 

13.610 

41 

18.300 

0.5128 

2388. 

50.88 

14.400 

42 

19.900 

0.4990 

2365. 

50.77 

15.600 

43 

21.500 

0.4912 

2345. 

50.68 

16.700 

44 

24.200 

0.4880 

2311. 

50.55 

18.400 

45 

26.900 

0.4874 

2281. 

50.45 

20.300 

46 

29.550 

0.4874 

2259. 

50.36 

22.100 

47 

32.200 

0.4874 

2237. 

50.25 

23.900 

48 

34.900 

0.4874 

2218. 

50.10 

25.700 

49 

37.600 

0.4874 

2200. 

49.95 

27.500 

50 

40.250 

0.4874 

2182. 

49.80 

29.300 

51 

42.900 

0.4874 

2165. 

49.65 

31.200 

52 

45.600 

0.4874 

2151. 

49.52 

32.900 

53 

48.300 

0.4874 

2137. 

49.39 

34.800 

54 

50.950 

0.4874 

2122. 

49.25 

36.600 


0 NSBL 

A BLAST, Ref. 3 
O Schlichting, Ref. 15 
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Figure 1. - Velocity profile for incompressible, laminar 
flow past a flat plate; angle of attack = 0® 
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Figure 4. 


- Heat-transfer-rate distribution for supersonic 
flow past a 12® cone, = 4.84; angle of 
attack = 0®. 
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(a) x/L = 0.10 

a 

Figure 6. - Temperature distribution across a laminar boundary 
layer, Orbiter Flight Design Trajectory; M = 9.49; 
angle of attack = 30.83®. 
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Figure 7. - DispjCwament-thickness distribution, Orbiter Flight 
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Figure 10. - Velocity distribution across a laminar boundary 

layer, Orbiter Flight Design Trajectory; M = 29.86; 
angle of attack = 41.4®. 
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Figure 13. - Heat-transfer-rate distribution, Orbiter Flight 
Design Trajectory; M = 29.86; angle of attack = 
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DEFINITION OF MATRIX ELEMENTS 


The individual elements of the coefficient matrices which 
appear in equations (28), (29) and (30) are as follows: 
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APPENDIX B 


The NSBL 
Progrva.m NSBL: 
Subroutine PIGYBAK 

Subroutine NUMIN: 

Subroutine DERIV: 
Subroutine FGES: 

Function TWOVAR: 

Subroutine EJOYCE: 


NSBL USERS GUIDE 
Program Contents 

code is made up of the following blocks; 

calculates conditions at the edge of th® 
boundary layer as well as s, 

called by NSBL - calculates the initial 
boundary layer profile using a Faulkner- Skan 
solution technique. 

called by PIGYBAK - solves ordinary differential 
equations using a multi-step method with a 
Runge-Kutta starter. 

called by NUMIN - calculates derivatives 
called by PIGYBAK - calculates first guesses for 
f (0) and g (0). 

called by FGES and EJOYCE - does linear double 
interpolation. 

called by NSBL - caxculates the downstream 
boundary layer profi.^es using finite 
difference methods. 
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Subroutine HOLIER: called by EJOYCE - calculates thermodynamic 

properties for air as a real gas. 

Subroutine DINT2: called by HOLIER - does double interpolation. 

Subroutine DINTl: called by DINT2 - does double irterpolation. 

Description of Input 

The definitions of the input parameters are given below. 

Figure B-1 presents a flowchart of the input sequence with the 

formats used. 

MM number of x-points to be calculated 

N number of y-points in the boundary layer 

KK equals zero for two-dimensional flow; equals one for 

axisyrmetric flow 

ISEN equals zero for variable entropy flow; equals one for 

isenrropic flow 

IPERF equals zero for real-gas properties; equals one for 
perfect-gas properties 

NTEMf number of elements in the transport property temperature 
vector 

NPRES number of elements in the transport property pressure 
vector 

DDSTR displacement thickness due to injection at the first 
station 


origin’al Page 
OF p(Mm quality 
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FLOW 

FINJ 

WMI 

WNS 

TEMP(M) 

PRES(M) 

VISSd.J) 

TCS(I,J) 

CPS(I,J) 

VISI(I,J) 

TCI(I,J) 

CPI(I,J) 

E 


value of the stream function at the wall at the first 
station 

value of the stream function of the injectant species at 
the first station 

molecular weight of the injectant species 
molecular weight of the stream species 
Mth temperature of the transport property temperature 
vector, °R (temperatures must be in increasing order) 

Mth pressure of the tiansport property pressiire vector, 
log^jj (Ibf/ft^), (pressures must be in increasing order) 
viscosity of the stream species at TEMP(I) and PRES(J), 
Ibf sec/ft ^ 

thermal conductivity of the stream species at TEMP(I) and 
?RES(J), Btu/ft sec 

specific heat of the stream species at TEMP(I) and 
PRES(J), Btu/slug ^R 

viscosity of the injectant species at TEMP(I) and 
PRES(J), Ibf sec/ft^ 

thermal conductivity of the injectant species of TEMP(I) 
and PRES(J), Btu/ft sec 

specific heat of the injectant species at TEMP(I) and 
PRES(J), Btu/slug 

convergence criterion on a boundary condition 
(typically 0.0005) 
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DELT 

EPS 

PTE(l) 

TTE(l) 

HTREI 

RGAS 

GK 

PRATIO 

TW 

ROVIN 

RDS 

X 

SEL 

SLE 

SELT 

UNIT 

ALFl 


step size in the y-direction (typically .05) 
convergence criterion when a variable step size is used 
total pressure behind the shock, Ibf/ft^ 
total temperature behind the shock, ®R 
reference value of the heat transfer coefficient, 

Btu/ft^ sec ®R 

gas constant of the stream species, ft Ibf/lbm °R 
ratio of the specific heats of the stream species 
ratio of the local edge pressure to PTE(l) 
wall temperature, ®R 

density of the injectant species, slug/ft^ 

radius of the equivalent body of revolution, ft 

streamwise wetted distance from the stagnation point, ft 

local entropy at the edge of the boundary layer divided by 

the stream species gas constant, S^/RGAS 

total entropy behind the shock minus the local entropy at 

the edge of the boundary layer, Btu/lbm °R 

total entropy behind the shock, Btu/lbm ®R 

equals zero if the initial profile is to be calculated 

in the code; equals one if the initial profile is to 

b*^ input 

initial value of the scale factor, a, at the x-location 
of the input profile 
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REXINT(l) 

S(l) 

BETA(l) 

r(i) 

THETA(I) 

CS(I) 


integrated value of the Reynolds number at the x-location 
of the input profile 

transformed x-coordinate at the x-locarion of the input 
profile 

pressure gradient parameter at the x- location of the 
input profile 

nondiraensional velocity (u/u^) at the point in the 
boundary Tayer 

nondiraensional temperature at the point in 

the boundary layer 

mass fraction of the stream species at the point in 
the boundary layer 


Description of Output 

The output for the NSBL code includes the input parameters, 
the free-stream parameters, the boundary layer profiles, and the 
boundary layer parameters- The output for the free-stream flow 
includes: 

M station number in the x-direction 

PE(M) static pressure at the edge of the boundary layer, Ibf/ft^ 
total pressure at the edge cf the boundary layer, Ibf/ft^ 


PTE(M) 
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TKETAE(M) temperature ratio of 'the boundary 

layer 

TTE(M) total temperature at the edge of the boundary layer, 
EMACH(M) Mach number at the edge of the boundary layer 
UE(M) velocity at the edge of the boundary layer, ft/sec 
URENOE(M) unit Reynolds n\imber at the edge of the boundary layer, 
1/ft 

REXINT(M) Reynolds number at the edge of the boundary layer inte- 
grated over the wetted distance from the stagnation point 
SEL(M) local entropy at the edge of the boundary layer divided 
by the stream species gas constant 
S ( M ) transformed x-c». ore i^nate 

ROVIN(M) density of the injectant species, slug/ft ^ 

BETA(M) pressure gradient parameter, 6 
FLCW(M) value of the stream function at the wall 

The output for the boundary layer profiles and parameters 
for M-1 x-stations Includes: 

M station number in the x-direction 

X(M) physical watted distance from the stagnation point, ft 

I station number in the y-direction 

YN(I) transformed n-*coordinate, n, at the point in the 
boundary layer 
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Y(I) 

ETA(I) 

F(I) 

THETA(I) 

CS(I) 

TAU 

CF 

HCOEF 

TREC 

QDOT 

STNO 

HRATIO 

DELS! 

THMOM 

DSTAR 


physical y-coordinate at the point in the boundary 
layer, ft 

transformed y-coordinate, n* at the point in the 
boundary layer 

nondimensional velocity at the point in the boundary 


layer, u(I)/u^(M) 


th 


nondimensional temperature at the I point in the 


boundary layer, T(I)/T^^(M) 


rth 


mass fraction of the stream species at the I point in 

the boundary layer 

local skin friction, Ibf/ft^ 

local skin-friction coefficient 

local heat-transfer coefficient 

recovery temperature , °R 

local convective heat- transfer rate, Btu/ft^ sec 
local Stanton number 


HCOEF/HTREF 

local displacement thickness, ft 
local momentum thickness, ft 

local displacement thickness with mass injection, ft 
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Figure B-1. - Flowchart of the input sequence and 
FORMATS for the NSBL code. 






WMI=WMS? 



Figure B-1. - Continued 











1 


Format (I10,4E12.5) 
IINIT,ALF1,REXINT(1),S(1),BETA(1) 



End of the input sequence 


Figure B-1. - Conclusion. 




APPENDIX C 


THERMODYNAMIC AND TRANSPORT PROPERTIES 


Calculation Procedures for a Perfect Gas 
Before one can solve equations (12), (13), and (14) for 
the independent variables F, 0, and C^^, it is necessary to 
evaluate the thermodynamic properties which appear in the coeffi- 
cients. The thermodynamic properties of the individual species are 
readily computed as a function of pressure and temperature. For 
cases of mass injection or when the free stream is assumed to be a 
perfect gas, the properties of the mixture are approximated with 
the ideal gas relations [10] once the relative concentrations of the 
stream species and of the injectant species are known. 


DENSITY: The density of the mixture is calculated using the 

relation 

p = (C-la) 

R T 

which can be written for a mixture of two gases as: 


p = 


« 1«2 


R T ”1 + ■ ” 1 ^ 


(C-lb) 


where, to have consistent units, the value of the universal gas 
constant, R, is 


1545 


ft Ibf 
lb mole °k 
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The transport properties (specific heat, viscosity, and 
thermal conductivity) for a given component are computed by linearly 
interpolating the input tables as a function of temperature (the 
transport properties are independent of pressure for a perfect gas). 


SPECIFIC HEAT: The specific heat of the mixture is computed using 


= C^, + C-(C^, - C ) 
P P2 1 pi p2 


(C-2) 


VISC0SI".Y: The viscosity of the mixture is computed using the 

relation [11]: 


M = 


h h 

^ °12 ^ Si 


where 


(C-3) 


S«2 
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THERMAL CONDUCTIVIT’. : The thermal conductivity of the mixture is 

computed using: 





where G. , and have been defined above. 

DIFFUSION COEFFICIENT: The binary diffusion coeffic-’’'‘nt is com- 

puted using the procedure outlined in [22] for diffusion in a 
multicomponent mixture: 


D. . 


ij F. F. 
1 3 


(C-5a) 


Where d is a reference coetixcient which can be approximated by: 

a. 5 


D = 


cT 




Further, 


F. 

1 


and 



0.»j61 



(C-5c) 


(C-5d) 


The consxant of equation (C-5b) was evaluated by fitting 
the equation to the tabulated values for binary diffusion coefficients. 

Specifically, the diffusion coefficients for mixtures of 
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AirrCO^, for temperatures from 520®R to 

12,000®R, as taoulated in [22], Using a value of c = 2,64 x lo*®, 
the calculated values were within 20% of the tabulated values over 
the entire range of pressure and temperature considered. Thus^ the 
equation used in the iiumerical code is: 


-j 


2.64 X 


•» A — ® 
..yJ 


TT 


“(■'TT 


(C-6) 


DIMENSIONLESS PARAMETERS; Tneo the dimensional thermodynamics of 

the mixture uave been c "ed, it is a simple matter to compute 

yc 

the Chaoman-Pubesin factor, , the Prandtl namber, — r-^ , 

e e 


pD 


and the Schmid c number. 
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Calculation Procedures for a Real Gas 
For cases where there is no mass injection, the free- 
stream species may, at the user’s option, be asstimed to behave as 
a real gas. When this is the case, the density of the gas is 
calcvlated by the subroutine MOLI^'R. The transport properties are 
computed by linearly interpolating the input tables as a function 
of temperature and pressure. 
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call -I. lien {hF(mT, PE rHT^],7ETN)#Z(H)#aEl6rNI«PH0E.(Nl,GA6MA(ftn A 1ii4 

UK ThTsSurT (^••(HTE Tht-mf T T *^Sft3l , (O A iftb 

E A f Ml asr;kT (Gamf A ( ~ * 3^^ 1 7u*RGAb* Tf (M ) •; fN) ) A Ifto 

kHACHC*l»tPE (NT/fcAfMT A lft7 

PlOGsaI OGIftCPE I^) ) A lft6 

VISE (Nl*T#*(JVAR fNPKES,NTFNP,PLnG# TE (n) #PRES, TEhr, viSS) A lft9 

THETA^f9)sTR(H)/TTF (1 T A lift 

ThETaE C'i)«TE(ft)/TTE ( n A lU 

UREnOF («)*WNOE ( 6)^UE (Hi/vISFT-) All^ 

6ft CONTINUE A 113 

GO TO tttft Alla 

9ft CUNT INOE A I IS 

no Iftft A lib 

SLM*^I«ft,ft A 117 

SELM>*ft,ft A IIH 

MlFfNlATTE (n*CPSfl)/3i,l79 A 119 

PTE(RT«PTEM) a l2ft 

TTEfM)*TTEf 1) A \i\ 

PECN)aPWATTO(N)*PTE(l T A l^^ 

E^aCnInJaSor: ( f (1 •/PRaTTUI^) (GK« 1,1 /GH )•!,)« (^«/(G6*l,n) A 1^3 







TE(«)aTTFCn/(l,♦^G*^•l,l•EHACH(n}•*^/2,) k \iU 

ET4fc(M)«TF(H)/nC(n A Ub 

TMfcTA^(M>«T-(M)/TTE(l > 4 12* 

HE(M)«TF (H)*CPS(n/i2.17tt 4 127 

PnOEC-^IbRE C“)/C7F (m1*S2« 174*WU4S) 4 12*1 

t4(H)«SQPT(r«K«32,t74*PG4S*TE(M)) 4 \?9 

UE (*i)«F4 rM)*FM4CMr«o k li^ 

PLOGS4LOK1 H (PE (••) ) 4 1 M 

VlSE(«i*ThliV4Q(spHFS#NTtiP,PLnC#TE(H) ,PPES»TEPP» VIS8) 4 112 

u4fcN0E(Ni -PhuE('^)*UP(«)/VISEC<) 4 111 

XiAiS CUnT1n>JE 4 lltt 

GU TU 16^1 4 11b 

continue 4 ll* 

PTEiaPTECn 4 117 

TTEI»TTf(i) 4 li« 

IE (iPtPP.EN.n GU Tf 1 Ik. 4 IIP 

C 4 lav) 

‘^fc4U a2-^. CSEL (►») ."**1 4 l**l 

C 4 142 

C 8fcL(^*> IS PE4f, IN. In IfcWHS OF S/M 4 I4l 

C 4 144 

C4LL mulIEM (HTE|,PTEl,2,TTEl»/TEl.STkl#MHnTEl,C4MTEn 4 I4b 

•in 12v^ HSl^MM A 146 

mTE(h)khTE1 4 147 

PE (M) =P»4 Tln(*^) •PTE 1 4 14ii 

CA'. L '^OLUM (mE(m>,Pe(M),1,TE(pT#Z(^ ) • SFL (h T . mhQE ( ^) ,G4mM4 (*• )) k lap 

t4('-)sSn»Trr,AMN4('->4l2.l74*MG4SaTE(^)*7fM)1 4 ibi^ 

L't(M)aSPPT( (HTE I -nf f a ) ) *s^*i62 . 6 ) 4 ibl 

E'*4CH(M>»0F (P)/E4(m) 4lS2 

PLUG2 4Ln(,UUPE (*^) ) 4 Ibl 

V lSk( ■^)rT#.nv4M(NPM4 b,NTk*lP,PLOG#Tk #P»ES, Tk«p, VlSS) • Iba 

UME NME (►•IsMMiip (^ ) *HE (P > / V I SE (M) 4 t SS 

call 4f*LTEM (HTE(f-),PTf (M),5,TTE('4),jTfM),SEL(‘^)»MHOTE(«).G4nTE 4 Ibb 

I (^n 4 lb7 

TMFT4M^TsT*»(M)/nF(M) 41 bM 

T-»ET4tf^TrTE(M)/TTF(^) 41 bP 

1 C jnT iNnfc k Ihk' 

(*n Ti' 4 I6l 

13 • ClU T Inhf 4 

C 4 1 6 1 

C SLfc(*’T IS WF40 \K n^TTS mF HTil/LP>^ M 4 164 

C 4 16b 

ak4U 42i4, ( Si,F (M) , MM> A ]6o 

ME4[» 42^, SELT 4 167 

C A 166 

<)U 141^ 4»1,mm 4 

SPL(4lsfSELT»SLE )# 7 7« , 1 6/404 5 4 I 74 

5LE(*^T*SLEC^:*2b^M ,4 A 1 7l 

1 4 < COnT I N(iE 4 I 72 

C 4 1 71 

C SLE(m) UNTTS 4Mfc IN FT LhF/SLUt M 4 174 

C SkL(^) IS TmE local SE/M A I 7b 

C 4 176 

hTfi.— PS(n*TTFl/i?,l74 A 177 

00 I , J M*| 4 176 

TTF(H)aTTEl A l79 

mTE ( ) bHTF I 4 16 ^ 

PE(M)«PkATinCMAH.i * 4 161 

rHEr4ErH)sEAP(l./rCP6M>*77S,l6)*(*SLk(M)^lZ»l7a*MGA$*ALUG(PMAr a 162 

1 ln(M))n jai 

TF(a)«THET4F(»')*Trtl 4 164 

PLOr,»AtrjGl^(PE('^)) A I6b 






VlSE(*4)sTMnVAti(NMRrs»^<1EMP«PtU6»TE(^)»PRES»TCMp,VlSS) A 164 

HE(M)«TEr‘<)*CP8n )/5J.!7a * 187 

EA(H)«SUPT(GK*329l74#RGAS*TE(^n A 16B 

THETA4(H)BT«i(M)/TTFl A 1§9 

UE(M)mS0RT ( (HTE t *he ( 4) ) *9lfii)6?,0) A 19P 

EHACH(M)sUf ( m>/EA(m) a t91 

PTE(m)sP((h)a (Cl •♦(&((• t.)/2»*EHACH(M)*EMACH(H))**(C4/(GH»Un) A 19^ 

HMOEC<)«(PWATiO(H)«PTEn/(Tr A 19S 

UPE4aE(H^4UHUE(M)*llE (H)/yiSEC**) a 194 

iSil CUNTlNUE A 19^ 

I6tf CONTINUE A 196 

print A 197 

PRINT 444 A 198 

PRINT 4S»S (htE(n),PWAT IC (“ 1 , T-« (n ) , tti)S ( H ) , K (N ) . ROW 1N( h 1 , -at , RM) A 199 

PRINT 284 A 244 

NSPACEAt'^T-R^^l/b A 241 

IF (NSPAC: .Ll,4) NSPAtE»4 A 2^2 

DO 174 isi.ttfSPACE A 24S 

print 29P A 248 

174 CONTInjE a 24S 

PRINT 484 A 246 

PRINT 474 A 247 

PRINT 484, f-,K(M),PEfN),PTE(N),THETAECHl,TlE(-)#tRACN(N) •N»l#NN) A 248 

PRINT 284 A 249 

NSPACfc»(62-**N)/6 A 214 

IF (NSPACf.LT.lM NSPACf«4 A 211 

DO IH4 I«1.nSPaCE a 212 

PRINT 29d a 215 

184 CONTINUE A 21« 

print 491 A 21S 

C A 216 

READ S54, NlT.ALFlrPf XlNT(n,S(l)#BETA(n A 217 

C A 218 

IF (IlNlT.FG.i) GO n» 194 A 219 

REilNTm*URfcNOE (l)«x( n A 2i?4 

S(1 )5RnoE( n*vtSE(l )*OEf l)*( (ROS(n*A8A)**2)*X(l)/(5.4A4RK) A 221 

194 Mat * 

print S44, N,^(M),OE(M),uRFN0E(N),Rf X1 nT(mi»SEL(**) A 223 

C A 229 

C OUrtNSTMEAM stations a 22S 

C A 226 

DO 244 ms2^nm a 227 

l)ELX(4)3i(Pi»A(N.n A 228 

RExlNTcM>»Hfc),iNT (N«i i^o.RtNHE CN«l)^iJRFNnt (MnA0F.LX(N)/2.4 A 229 

PRINT S44» *^#X(4),uF('-l,ilREN0E(N),Rf X4 NT(n).SEL(N) A 234 

OnSTR(Nis(ROVTN(H)*(prsrN)**f(K)/(RHOE(M)*uE (N) )^RuvlN(M*i)*(Rna A 251 

1 (N«n**A«0/(RH(jE(N.n*i|^(H*l )))*UEtX(H)/(2.«)A(ROS(N)«*KK)WDaST A 2>2 

2 R(N.i) A 253 

C A 2 59 

C EVAtUATION OF S A ?3b 

C A 236 

S(N)aS(H»1 )45i.b*(RHOE(N«n*UE(N»l)AVlSF(N-n*C (R0S(Nwl)A*88l#*2 A 237 

I )4RHl)E(M>*ut (91 AVtSE {N)*( CRoS(N)**8*i)**2n*0eLA(N) A 238 

240 continue a 239 

DS(naS(2)*SCn A 294 

IF (IlNlT.Fn.n GO TO 214 A 291 

8ETA(nac2,ti*rS(2)A$(l))/(uE(2)AUE(nn*i fuE(2)«uE(in/(S(2)»S(l)) A ?42 

1) A 293 

214 RETAisHETAf 1 l A 299 

00 224 Nsp,MMAk A 29S 

0S(Nl«(S(N»i)*S(N»m/2,4 A 296 

8ETA(Ml«(2,4«S(*i)A(UE(ri4 1)«IJF («•! )) ) / (Of ( n) * ($ % ] )«s (N*l ) 1) A 29 7 
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ii<» CONTINUE « 2«« 

C « 2«« 

c ev«li)«(.tiqn or Fcc *t the hall * in 

c * 2S1 

OU 2Si< Ha2,MH A 2^2 

f INJ(H)ar|Nj(M«nA(0OVIN(H«| T*(NOS(H»n«*«tOAHOvlN(H)*(H|)S(H)A* A Hi 

t KO )*I)EIX(H)/2,M A 2S<t 

FLC«i(nT«»FINJ(H)/(SO0T(?,*(*»CH)) ) A 2S5 

2J4 CONTINUE A 2^6 

C A 2b7 

PRINT 281^ A 2S« 

NSPaCC*(N2*HM)/6 a H9 

If (NSPACE.LT.a) NSPacEpp N 2PP 

UU 2<»H I*I«nSPacE a 2N1 

PRINT 2RU * 2P2 

2ttt« continue a 2«1 

PRINT SH A 2»P 

PRINT S2P» (N.AlNT.SIMl.RilWINlHT.eETAIHI.FLCHlHT.Hml ,HMAK1 A 26S 

C A 266 

C CALCULATING TmF FLU" PARANf TFRS AT THE INITIAL 8TATIIIN A 267 

C A 268 

IF (IINIT.EU,^) GO T». 2S«I A 266 

*1-F«ALFI a 270 

C A 271 

KkAO ia.i, (F(n,THETA(n,CSm,ls|,N) A 272 

C A 27J 

tio TO 260 A 27u 

2S0 '<61 A 27i 

oElsoE(l) A 276 

HTFISHTein A 277 

T"lsTn(n A 278 

TTlaTTEd) A 276 

CALL PICTbAK THTFl) A 280 

26V» call tJUTCF. A 261 

c A 282 

270 FiiRHAT (71101 A 2«i 

280 FORhaT (IHI) A 28U 

2*0 format (/) A 2«S 

38R format (0E12.S) A 2"6 

Six format (^>ei2,5> a 267 

J28 format (b6i,2a(lH*l,/,86A»2«H* TmE InPuI PaRAmFTFRS *./,S6X,26(1ha A 266 

I),//) A 266 

53U format (a6X»9MTHFRE AHk,I3«26M STATIONS IN ToF * UIRt CT ION,/, «8*,R A 2RX 

IHTHMt AwF,I3,26h stations Tn Thf y OIRFCTTOn,/) a 29i 

360 format (SX,10nH*),ll6MY')0 HAVE INPUT AN ERRUnEUUS A", A" MUST EOU A 2R2 

Ml 0 For a TmO'DImEnSIONAL CUnF IGURATIUN, or 1 for a 800V of *REV0 a 2R3 

2LUTI0N-) A 2R6 

3S'< format (S6X,«ShThE 6(>0V is a THOaOlMENSlONAL CONFIGURATION) A 2RS 

36<A format (S2X,32hThE HoOV IS A HUOV OF R( VOLUTION) A 2R6 

37R format (/,37i,OOMThf MOlFCuLAR nElGHT OF THE INJECTANT IS,F6.2,I3H A 2R7 

1 L»P/LHM-MHLE,/,3»*»iMMTHF molFCULAR «EIUHT of TME SPECIES IS,F6,2 A 298 

2,1 3P LBm/lNm.holE,/) a 2RR 

380 format (b7x,2lHTHE GAS IS A REAL bAS,/) A 3«lvl 

390 format (S%x,26hThE gas is a perfect gas,/) a 3«>1 

68R format (62A,9H00STR(na,Fl?,S,18M, FLC" ( 1 )a,El 2.S, IRH, FlNj(t)a,Et A 3R2 

12, S,/) A 3R3 

618 format (63x,2hE6,E12,5,7m, 0F.LTs,E12,b,6H, EPS6,El2,b,/> A 386 

628 format (6F12.ST A 3^6 

63R format (tX,|l6(lM*)) A 186 

668 FORMAT (|HR,17x,6HMT£r"),lAX,9MPRAT10(M),11X,bHTH(M),|«X,6HROS(Ml, A 187 

n3«»«Mx(M),|sx,8HRoviNrMn a jmh 

6S6 format (16X,El2,b,b(3x,iH*,3x,E12,Sn A 589 
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fUt»M*T (S8*, n^'I'OTP'lT n*T4,///) 4 JH 

uTk' FUPM4T (2lx, iHXf 1 1 (fSMPMM) ,8 i,4mPTF(M) ,7i|9HrHET4t (M) (bitb 4 ill 

lNTTE(M),8i«,enfcP*C- n * 5>i* 

MSP PU0M4T (t7x,IS,6Elu.5) * MS 

ttflp F0BM4T (29», JMM,T**PX*(H),lSX,5MI»fcfN),TX»9MU«fcH0t(*<)*5«.PHkE*lMTfM 4 3J9 

1 )>SX, 6 HSELfnn S 3 )^ 

SPP FU 0 H 4 T ( 26 X| IStSElUtS) A Mo 

SIP FU 0 M 4 T (, 2 PXt H* 4 «*A,aHX (M) , tPAfanStM) f 9 X«aHH 0 vis(P) f 6 X( 7 HaET 4 (M ) , 7 4 3 l 7 

IX, 7 HFtC*<(H),/J 4 il« 

S2P FHPM4T (2bX,lS,SElP,^) 4 3)9 

S3P FUPM 4 T (I)P,mE12,S) 4 un 

sap F0HM4T (SFIP.P) 4 32) 

C 4 Hi 

tf 4 323 



SUUHUUTlNC PIGVS«K (HTEI) 

C 

c suhhoutine pigvhak Calculates the initial phufile 

C USING The HLAST PHUCEOURE, 

c 

niHENSIUN F(S0)( GISPIi FN(9H)« GNTSlAIt TN(SP)» ETAO(SP)f ETA($i4), 
I THTTA(5«) 

OIHENSION V(S)f FPA(S), GAf}) 

OIMENSION OV(S), PIS), TEIS), rPN(S.P), DEP(S) 

DIHENSION EHIS), OPSAVE(S), FPH(?«>iA) 

EXTERNAL UERIV 

COHHUN /NE»/ F,TMETAISm),ALF,NV,KK,MM,C8fSP)i»t6IA).IPEHF#NUHPl 
COMMON /CONST,' Al#A2,Al,A«,E,EFSfOELT, *'’C' ' 

common /main/ Tmi,UEI,ITI 

c 

C CAUTION *PNbH,7* initial profile USES PeK. eCT gas RELATIONS 
C 

ItPMAM 
AE T«n 

AI«UEI*UEI/(a.HAMlEI*J2,l T«*77»,) 

A2tl,iA>Al 

PRA/,7 

A3«I,0/PW 

AaauEl*UEI*(PH*t ,P)/( hTEI*PR* 32, t , ,) 

PHINT At,A2,A3,AU 

AAtri.ll 

ODaT*I/TTl 

CALL EGCS I00,CC,EE) 

Ci C T n 2 r > 

1.1 CC»,999*CC 

PmInT 21/1, A*,^b#CC#l)l)»t^ 

INaN^'*! 

V(l)aAA 

r(2)88H 

y( J)»cc 

V(u)at)f) 

y(5)aEE 

fiPSAVt«yc5) 

FPPSAV«y(5) 

FPPk'iavc n 
GPil|aV(S) 

FPPOSavn) 

i;pos*y(5) 

ilal 

NaS 

an CAtL NUHiN (N,D£LNT#t)V,P,Tg,EB,yPP,y,EPS,r'PSAvt,KEy,aE»lv#n 

T.n. 

NmS 

laj 

CAUL NoniW (N,DELT»T,f)V,P, TE»e«#yP»#y,fPS,nPSAVt»KEY,OfRlV,2) 
in CALL numin (M,DELTf T>0 v«P, rrfEP,yPR«y,(P6#OPSAvefKFyfOERlV«3) 

IF (lERW.EQ.n 60 TO in 
IF (AHS(Y(3) i.lT.E) 60 TO 
lalf} 

GO TO Sn 

6n IF (Ad$(l,nnnn#Y(2)),GT,F) GO TO 70 
If (AHS(i«nnnn»Y(a)),LT.E) Go to inn 
70 FPA(inaY(2) 


H 1 

H i 

H i 

R a 

H i 

B 6 

B 7 

B B 

H 9 

H in 

H U 

B \i 

B 13 

H la 

15 

H Ih 

H 17 

B IB 

H 19 

H 20 

B 21 

B 22 

H 23 

« 24 

H 25 

H 20 

H 27 

B 2« 

h 29 

B 30 

B 51 

H 32 

B ii 

P 34 

H 55 

b 56 

H 57 

H 5B 

B 59 

H 40 

H 41 

R 42 

b 43 

H ua 

H 45 

6 46 

B 4? 

B 46 

B 49 

B 56 

R 51 

B 52 

b 53 

H 54 

R 55 

B 56 

R 57 

H 56 

H 59 

H 66 

b 61 
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a 


[ 


r. 
[ 
r 
[ 
r 

i 


UA(ll)«v(a) 

lialUl 

If (n,EU«2) GO TU 

IF (II,NE,3) GO TO 

Y(1)BAA 

V(2)bOH 

Y(3I«FPP5aV 

Y(4)aD0 

Y (5)B1 ,^i95*GFSAVE 
FPP03«Y(3) 

GP^3«Y(«i) 

FPPOS«Y(5) 

GPOSsY(S) 

GU tn au 
614 V(1)«A* 

Y(2)a6H 

V(3)«1,«90S*FPPSAV 
Y(tt)sOD 
Y(S)sGPSaVE 
FPPw2*y C5) 

6P02«Y(5) 

FPPUSsY(3) 

GPOS»V(S) 

GO TO 40 

90 UFPABFPAM )«FPA(2) 

0GABGA(n«GA(2) 

OFPAiaFPA(t)«FPAr3) 

OGAlBGA(l).GA(i) 

l)FPP0BFPP0)»FPP02 

OGP0aGP01*UP03 

A6*UFPA/OFPP0 

M|«DFPAi/nGP0 

ClBt,0000«FPA(l) 

aiaOGA/D^PPfc^. 

F i«ocAi/r)r>P0 
Fia|,0000.GA(n 
OF.FPPa(Cl*fcl»Ft*Hn/( A6*F 
i>EGP«(A6*Fl«Dl*Cl^^tA0*e J«oi*h1 ) 

FPP«FPP01*OEPPP 
GpaGP0ltOEGP 
Y(l)aAA 
Yr2)BBri 
Yl5)aFPP 
Y(4)aoo 
y (5)BGP 
GU TU 30 
l»0 LJnTInuE 
YCDbaA 
y ^2)BBR 
y())BFPPO$ 

V ( a)«ou 

y(S)BGPos 

CALL NUHIN (N,r)FLT»T,nv,P,TE,fcH,yPW,y,EP5,OPSAVt»Kfey#UE«IV,n 

Tb0, 

was 

ULSTia(y (4)«Al*(y (2)«*2n/A2«y (2) 

THTAl«(y(2))*(t.0«y(2)) 

THfEiacy (2))*o ^0«(Y(2)**2) ) 

1 1 ia| 

FP«Clinayf2) 

GUll )by(4I) 
feTAOdinaT 


B 42 
H 63 
B 64 
H 6b 
6 66 
H 67 
B 6B 
H 66 
b 70 
B 71 
B 72 
6 73 

6 74 

8 7b 
H 7a 
H 77 
B 76 
6 79 

6 60 
H 61 
B 62 
b 63 
P 64 
P 6S 
B 66 
6 87 

H 66 
H 69 
H 90 
b 91 
J 92 
B 93 
b 94 
B 9S 
B 96 
H 97 
H 96 
6 99 

b 100 
B Ml 
B 102 
H 103 

8 104 

H 10b 
6 106 
B 107 
B M6 
H 109 
B 110 
H 111 
H 112 
R 113 
B 114 
B 11b 
6 116 
B 117 
B 116 
B 119 
6 120 
H 121 
L 122 
B 123 


* . 




CALL NUrtlN C»^#OtLT#T^OV,P,Te,EW,''Pk# V.t P8#rPSAVE,^EY,ueHly,2) 
ou8r2B(r(4}*AM(y(?i**?))/A2-r(2) 

THrA^BrY(^))*rl,«9•y(2)^ 

THTt2»: M^))•(l•^•cy(^)••^i) 

AUELSTai^,8*0LSTt ♦OLST^)•UELT 
AtHf TA«(»1,5 *ThT*1^TMTa?)* 0ELT 
ATHETe«{lfl,5#THTU4TMTE2)*0ELT 

Ilk) CALL NU IN (N.,[>ELlf T#0VfP#Tr#£NfypK#y,EPS»nP8AVE>KEy,0ENIy»3) 
DDeLST»c(y(«).AU(yc2)**an/A^-yc2) )tpti i 
i>THETA«( :yr?n*(u^*y(i:)))«OLLT 

i)Th|tc»c f y ( 2 ) )•(! •*d-(y c2^ **2) ) )*del r 

A1)EL$T«ADLLST4DI)ELST 
aTmetasaTheTa^othE fA 
ATMfc TE«A fht TE ♦DThETE 
Ill>IlUt 
ETADvII l)«T 

Fp^ain«y(?) 

&nn )sy(4) 

If CNCOUNf.Nt.i^) CO Tu lik) 

If (FP*-ain-M.Q9) I 

CALCl»LATr)N Of AlP 

I2H ALF»*AL0G( 1 ,/FLOaT(In, )/F’ 4f)f 1 1 ) ) 
nCOUnT» 1 
15^ CuNTlNiit 

IF (AHS(y(jn,LT,E) GO Tf‘ la.) 
bi* To Uvi 
14^ CO»^TInuE 

STk.A** STATION TNaNSFU^^m •IN* NO, UF SlATIONS 

00 I5k» I»l»Nr 

yMnsFLOATfl-l VFLOA (IN) 

ISO COnTIniJE 
J»1 

t Ta ( jl«ALnf,( 1 ,k»*YN( J) n/ALF 

ou I VO isum 

!f>r^ IF (£Ta( JI-FTaOU ) ) i7rt,iv/,19yi 

17^ 1I»I*1 

AFS(fcTA f JI-FTaDU m/(ETAO( I)»£TADf in ) 

Mj)sFPfcfnn^(FP«rn»FPw(in )*af 
gn(j)sg( in^(r,n)-c(M))AAF 
CS(J)*1,0 

THtr4v(Jl«GN(J>»Al*F( !)*F(J1 
JaJtl 

IF (J.bl.lN) Ou Tu 180 
fcTA(J)BALf»(,n,0/n .O.YNf J)) )/ALF 
bO T ) 160 

180 fcTA(J)»l ,0E 10 

190 COnTInije 

F(IN^n«l ,0 
rN(iNtn*&i n n 
Theta , iN^l )si .o^At 

CSUNtnal.O 
HE TURN 

^00 F OHM AT (^K,: 'lAtB#F9,S,2x»3HA28^F4,5#2X,iHA*«,F9,5«2X« Ar«AttS,FQ,S) 
210 FOPhaT (2x,3haaB»EI?.S»2a»3hh- i.»fc l2.S,?y , inu.%. El2,H,^x,3»^O0»»fcl2, 
1‘>,2X,3hee*#E12,S) 

En:) 
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H 
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H 
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1 40 
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1 44 
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1 49 
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ISH 

H 
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juHMOUTlNt »U(S (DU.CCftE) C 1 

/nuS/ HfT*l C i 

OlMfNSIUS KETiS(S)« FI»P(S,«), CP(S«<0 C J 

o*Ta G/.24..i«.«.,oS/ C 4 

.'*r« i»ET4S/.iiIS«.I,.2>.3>.4/ C % 

04T4 FPP/, 4179121, .«74S6T^,.STSli«7,.6(iaS224.. 7451, ,«lbbr91,. 49243 C • 

147.. 5941555..45178 39,.76.->P43,.4S«4/(I3..515447 3,,617««69,.T279P41, C 7 

2. 7494819. . 44B178P,, 5254897, ,6274183, .7174249, ,7989913/ C 4 

04 T 4 CP/. 2419834,. 2474471, ,254>»55i,, 2424474, ,248248, ,2335142, .2385 C 9 

1798.. 24444P9,.?528482..2581825,.2121744,.21449.15..223541,,2287342. C 14 

2. 2331344. . 1988147. .2428812. .249193.. 2141324. .21811;27/ C 11 

*■8*5 C 12 

8(i«4 C 13 

CC4T*OV44(Sri.NG,RET81,l)U,8F14S,G,FPP) C 14 

EEsT.OV«4(48,K5,hCT*1 ,0(;,H3 rAS,U,i,P) C 15 

“EUiW't C 14 

C C 17 

bs. C 10 





SUHkDuriNE KUHIn (N,orLT»T.nVtP,TF«C(7«VPH,r)CP«EP8tOPSAVEtKEv,OCmV 
IflMDX) 

COMHON /const/ At«*?»AS,Aa,e#E^*18f’N 
C 

C ••— •C€Nf»4t NUNENIC4L 1NTE6«4TIDN «ITH SInGlE«STEP EB«0« AnaLYSI 

c 

01NEN8IUN a(4)« R(S), Ov(n1, P(N), TE(N)» EM(N) 
dimension rPPfN,4)# UEP(N), OP6AwE(N) 

C 

C ••••«• 

c 

c 

EXTt«N4t UE«IV 
C 

C <.•«•••••«•• 0AT4 

c 

6U TO iNhk 

SiM«2a.P 

/EPOStf«/ 

Am»«p,«d/six 

A(e)«37.-^/sxx 

4(i)*«S9,fc“/SIX 

A14)»SS,P/SI* 

b(l)sui4/six 

W(^)*•S,M/SIX 

b( 3)*19,^/sn 
cMU)s«4(n 

Si 

Tno»^,4< 

Ni»a 

M5a^ 

«ETUPn 

c 

c ••••*••• SET UP INTEGRATION 
C ••••«••» ENTwt NfSET 

C 

assign Tn JPLS 

KOPst5 

>^uunT«o 

OElHV6»0ElT/SIx 
OEL bY^»UEuT/Ti*0 
0(5 $ta J»J,N 

UPSA^E(JT»0EPTJT 

t)V(J)sDEP(J) 

5>A Continue 

CALL OERIV (T,nv,P,6fN) 

IF (TEPK»E^J.n return 
DO 4^ J«1,N 

yPR( jf Nn*p( j) 

ER(J)«/ERO 

«*» continue 
return 
c 

c ••••«.«• entry-* >OInT Enw ONE MJRErICaL INTEGRATION STEP 

C 

SR V»T 

T»V4DELT 

»<OP»J'>-TP^I 
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s 

4 

s 

p 

7 

tt 

9 

IP 

U 

li 

IS 

14 

15 
1 o 
17 
IS 
19 

dl 

^5 

^4 

PS 

po 

P/ 

Pb 

P9 

SR 

SI 

SP 

ss 

54 

iS 

iP 

57 

Sb 

59 

4R 

41 

4P 

45 

44 

45 

4ft 
4 7 
44 
49 
SR 
SI 
SP 
Si 
S4 

ss 

56 

57 
Sft 
S9 
6R 
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o n n 


Go TO 1 Pl5 

•••••• «U^CE-*<UTTA PkOCEOuRf* 

vvsv^oELava 
OO J«1,N 

OV(J)«OEP(J)^YP»CJ#Mn*oEL«vi 
74* Continue 

Call OERIY (vv,ow,p,3,N) 

IE aE««%E«,n «ETuRf- 
OO JPl#N 

DV(J)«DEPf J)^P(J)*OELHt? 

H0 Continue 

Call OEwtv (vv.Ovf TE#<i»M 
IF (lEPR.EJ.l) weTUHs 
Ou 9^ J*t»N 

Ov{J)«OEPU)aTE(J)<»OELT 
Tt(J)«2.W*CTf (JUP(J)I 

9M Continue 

CALL oe^lv (T,Ov,P,5»N) 

IF (lEWH.EO.II met JWN 

OU J«1#N 

OEP(JT«ntP(j)^ntLHVb*(Pf j)fTE(JlAVP»(J.*^t)) 
OY( JT«0EP(JT 
i*ui continue 

call OEMIv (T,Uv,P,b#N) 

IF (lEPM.eo.n nETJMn 
OQ ltk< I*UN 

IIP continue 

AOliNT*KOUNT4l 


C 

C •••»•• CMfC*< the nu**HEM of InTEGMATION STEPS hadE «V m«k 
C ►'ttriCEDUME 


IF (nCUnT,lT. 3I r,(. If 
assign TO IPl5 

HMsHa 
h4sm3 

IF (‘<ET.EL,P) METumn 
IF (‘'HP.le.S) meTuwn 
^lNC«0 

OU J«l»N 

IF (E«(J),GT.tPS) GO ro iSvl 
IF (£M(J)*14 ^i^v4vi,4I,lT,EPS) KINC»KInC*1 
13^ continue 

IF (KiNC.tQ.N) GO TO 19« 

OU JS|«N 

OPSAVE(j)soEP(J) 

CONTINUE 
MEruPN 
ISO CONTINUE 

IF (FOP^EQ.a) Gh to 164* 

TsT^DElT 

Ipo 0ELT«0.S#UELT 
OO 170 J«l,N 

DEPfJ)»OPS*vt U) 

170 continue 
GO TO ?.A 
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SUrtWOUTlNf ptWIV ( t , Y ,!•< ' 

rU'^fcMSi'iN vrs), nY(*>) f '' 

Ci>^Mi)N /DOS/ HfT*l ► -> 

common /Const/ * 1 • , A 5, au , ^ , t>^S • l T , U m ♦• 4 

X*( (Y( 4 )»( Al^Y (^lA*?n/ A/> t ^ 

TF (x.tE.i^,) no To K »■ O 

Al3(X*#(*i,»,/Sn t / 

X^a ) * r »** ( - l , Y?S) W ( ( V (SI / A / ) •^, * ( A I / A? 1 * Y I / ) * V r j ) ' ^ f. 

0Y(n*v(^) 

nv(^)»v(i) ^ 1 .’ 

f)Y(3>a-n,/xn*(fV(iw>^)#Y(n^^tT4i*(«i.y(^*Y(/in ^ ii 

l)Y(tt)SY(S> f 5^ 

i)Y(S)««(i,/fA3*xt))*((Y(n4A5*«/)*T{'>WC(xi*v(i)«*/l*tin*Y(/)*y(.i t 1 i 

n)*Aa) ^ I- 

k^FTilSN ^ i3 

1.1 ^ 

WetUKM ^ 

c ^ t- 

fcNO '' * " 
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• ^ 

• i 
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iUHWOUTiNr FJnvtF 


SUmMOUTmF FJijYCF OSfcS MNtT^ 

To C*LC»JLATF ThF OO^NSTwFAn PN<lFTLfcS 

OiMfcMSiUN ThMaKSi:), CSUb^T, OfLri)S(S^). onrsu 

Tf^V»), FF(Sfc4)# Fi.Ct'JoU FLClCSfcl), T hf T ( S/. ) , CS^Cj.n. FLC^fS 

DlHtNSlO'i CMi(Sk')# CT*SfS%M, T# vISSfSM, Tcnx|(S^i), TCOv*SCS 

1“)» OirFC(^«T)f kKS*M. 

TCOSM(^p) ^ UnOm (Si? ) * 0**^ f • T • .>L rs ’ ) , ) 

• 'l^t'^SION FTAfikl)* OFTA(*'»n# V(S.)» iJfctbTfS»>)# I h«^i ( «>^ ) , i ( ^,M 
tJl ife'^i>ION OsTAKCS^)# orssTKfS^T 

I' '*iSIU*« A (SO, it 3) ^ ^ f S^'f i t # C (*'«’, S t (Sf, )t 1 ) t *<( S.' t i# ^ 1 • 

lHAH(SOtiti)t t^Ahl^^'^T^'-^flfiTt A»MrAH(S/,i,|), A ^ (Si, 3,1), o(S*^,itl) 
if i)AG(so,itn« i»(s*'*ti,i )t H-i(s<,^,n. -‘(soti.n 
I) l TAifcfSOlt UOI- T t f ) , ST^ <TC^S*'), j j* T 

I«(S15) 

>lMf Tc^**‘(5vi)t nN,»T(S.M 

dimension «K^) 

ro-iMON /Pkj)P/ TfeMP(a;d),pk|S(a),M£Mp,sFN(>s 

CDhmON /u^^/ F (SO)t r A(S/),ALF,^,IO•,-M,CS<SOT,»,IP^.M^ 

Cu*<kliN TTF (to^T t»»^T # -^•^St»'0S(60) , OOSTP 

CO*^HON /nF m %/ FtC"(0O>, TMf. T At ,klJVlN(801 ,hf TA(60) ,«JS(F»»^) t ‘T , 

IU£ (60) tkKUF (601 

/MF«a/ PF (60) , Tmf T A-.(6* ) , TF ( e». 1 . V ISr (6ti) ,HT«f ► 

CljOHUN /NFwl O/ KF. Vt V6C (1 6<r1) , V 1C ( t 6i' ) » U S( ( 1 6 .' ) , I C i f ( 1 6 .* 1 , L ^ i C ( t 6 ‘ > 
i»C*»5r. (i6i) 

•JtLTsl, 1/FtUAT(N*t ) 
half 30 
kONl«/stS<iS, 

Fl3(i«^T/^6,»t)^«t:i6l 

FSs(**«S/*?6,/)#*,a6l 

MSI 

1 

t)0 10 ISitN 

Y>4( 1 1SFLHAT( I-l )*KFL T 
F I ( 1)SF n ) 
iMt TAI ( I isTmF TA( I 1 
CSl ( I )«C6( 1 1 
10 ConTInUF 

F2( 1 1«F I M > 

F2(i?)=F 1(21 
'>0 20 Xsi,N 

nFLCOS( n=o,n 
UFlCSI ( I )»DFLCOSn 1 
TUlsTTfcf^lATHtTAd) 

2^ (.(f.^iNue 

ow 30 I^l,^J^4 

FF(nsF(n/(ALF*n t r i-n *i.u 1 1 1 

30 Cow r I SHF 

tV4L'»ATl,)N l)F the STMpAM Fit^r T Trr.,F lC 1 

FLC n )*FLf.**(H) 

^lCC<2)=FLCnW(F(l)/ALF^P(2)/(ALF*fttO-l'tLTm*OELT/2,* 

C('j ao 1<3#NN#2 

FLC ( I )»FlC n*2)A(FF ( I.2 W**,o*fp ( I«i ( n 1 ♦htLT/ 3,0 

IF ( ( I ♦! ) ,nT,Nk) liii Tn yo 

FLC(iAn»FLC(T«n^(FFfi*n*4,OAF»^nwFFn4in*i)FLT/3',o 
ao CumTInuF 




DO SD l»USN ^ 

FLCtcn«riccn ^ 

COMTiNUP ^ tea 

N5PACt»CSa-N)/«> If 6S 

IF (NSl>4Ct«lT,^) MSPACr*/ F »6 

PUINT 69«t f 67 

00 6D IsunSPaCE F 6H 

P»INT F 69 

61^ continue f 7v 

pmint 71M F n 

PPINT 7^V» F 7^ 

P«INT 75H, (l«F(l)»TMETA(n,CS(|),lsl,M) F ^5 

l«l F 7a 

PtO&*4tU(;ifc'CPF'"')1 F 7b 

CPS(I)*T-»OVAWfNPMFb,f«TtMP,PL(’t»* tn>#PI»FS,TF«M,rPbr) F 76 

vi5S(n«i»uv4w(NP»Es#*-re^‘P,Pi ua,T(n,pi-i$*Tfc^p,vsc) f 77 

TCONSC n*7i*»0WA»CjP»fcb#'^TF«P»PL0r;,f (1 ),PPtS#TFMP,lf SC) »• /»* 

IF (KbV.(:(I.n GO to Til F 79 

CPI C I )sTi»0VAW(NPMES,9tt'^P#FL0(i, I Cl)»P»<FSf 1FHP,CPIC) F 

UhnI ( naT**0W4«(^P9Fb,Nlt*^P»PLnK.T( I ),PF»FS,Tt«P#TriC) F rtl 

'^ISKl1aT»*nV4P(NP»tS#NTF^PtPUift# TU >,PPtb,TtHP, vIC) F 

Gn TO ftp F fti 

7rf CPIU)aCPS(n F fta 

TC')Nl(n*Tcn‘^S(n F 6b 

visirn«viss(n f 6^ 

h** If ^Pt^^F,£Q,n CO TO 9P F *»; 

CALL (H,PF(ft)»T»TCn#/»FM#PrtOSf n,C) F Hh 

9H0^( I )sOmOS7 I ) F 09 

U|> To F 9/ 

9P 9hOH(| )SPF v*|p* I 7a* T r 1 )*(*»>^S*CSn ) • ( I • ) > ) » 9i 

lOP continue F 9<» 

cpp(i)scpi M)4cscn*frps(n*rpi(i n i 9> 

IF (C^C|)*l,U) F 

IIP *b(n«fCb(n*p«n/(Nf*sFCsm*(»^fti«i^'<s)) f 9:, 

*i(n«i.*i**srn f 96 

If (1*1(1), F0,P,P) go to F 97 

ui^n)3((i,p*so«T(vissn )/visif n )•( (^^i/«*-3)**p,^b) )**^)/(^,ft^ft<l* f 

ISOPT (1 •♦4MS/6NI) 1 F 99 

U^l ( 1)S(( 1 .P^SOkt f VISK 1 )/MSS( 1 ) )•( (»«»’S/»M^♦*4,^b) F l/P 

1 SU><T n *F-lNl/**MS) ) F 

vis«( naviss( n/r 1 ,o*ui<?r n*<u i)/*bcm*vibiM)/( i (n*i*b( i ) f 

1/Xl(t)) F 

TCONH( l)sTC(D4S(n/( UP*1 .^hb*Gl^( I )♦*! n )/«S( n )♦TCuNI ( I )/(l,P*l ,p F 1^4 

iPb^G^l (D^FSCn/xId )) F l>^b 

un TO XJP F 1/6 

12P V]SN(n3VlSS(l) f 1/7 

TCON^d )sTCU9S(n F 1/6 

liP continue F 1/9 

TAUSPH 06 ( n*VXSH( n*UF fM)« .E(M>*r90S(Nl*ftKH)*ALF*CF(^).F n ))/(0ELT F 1 U. 

l*(f2.P*S(H))**p,S)) F III 

CF 5t Aii/(P,S*9MnF (N)*OF (H>*I.IE ) F 

HE AT*«TCONMf| )*rri(H)*oE fH>*kHnM(| )*7 hi5(h>»*6h ) *ALF ♦ ( T«fc T a F Jli 

^A(l))/(DfLt*((^,p*s(ft))♦*p,sn F U9 

90(1T*«HEAT F lib 

I«N F ll6 

PLUrt«ALUGlP(Pt(ft)) F 117 

CPSa>»T60VA9(NPRfS,NTEftP#PLUC#Tn ),P9FSf TFHP^CPSC) F H» 

VISSC I )*TNOVAH(NP»tS»OTtwPfPLOr*, T(1 ) ,PWE S, Ttf^P # VSC) F 119 

TCONS( I )»TwOVAP(NP9ES#Nlf MP^PLnC^t ( 1 ) ^ PPfc S • TE«f # TCSC ) F l^W 

TCONft(I)*tCOMS(n F l^l 

VISH(I)»VISS(I) F li>^ 

CP'^(I)»CP5(1 ) F l^i 
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0 

P 


pwNot n»visM(t)*c»>^(n/TCoMMf n p 

MtCaSQRT fPrtNn(^) ) f’ 

T»EC*HEC*tTTF f ) -Tfe (Mn ♦TF { m) f 

S!NO«tiOO?/(»HnEC»*)*Wf. fH)*CPM(Nl#c TWtC*fMfc T A^(->*T IF f*^) n F 

PRINT 7a«», T*u#cr 9’ 

print 7514, TREC#UD0T,STNP F 

F 

ULOOP*^^ F 

14^ COnT Inm^ F 

15«4 CONTINUE F 

4LnOP»NL'UjP^l F 

IF (rU^OP.KT.ISU) UU to F 

PtOG»4LOGl^rPE(>^^ J F 

DU t«I»N F 

C F 

C EVAUUTIHN of TNf Thermo PRuPFNTlt> t»F LU'iPu>»t >*T S F 

C F 

CPS(n«T«0VARfNPFtE5,NTE'^P,PLtU., U 1) , PRt S . T L ^ P . CP >«n F 

VTSbC n»T^OWAR(NPRF^," TF‘lP,PL‘iU, T ( p , PnF S » T 5: • v F 

TCONS<P«TiiOV4R(NPRES,NTEHp,ptijr,, T ( 1) , PPt S , T h^P , T f Sf ) F 

IF CKET.fQ.n CO Tu IfeP F 

cPi(nxTwnvAwrNPMFS#NTE*^p»pi Mu,T(n,PRFs,tF^p,CPio f 

TCUNUl)*T«DWiR(NPPFS.NTk»^P,PLUf;,Tm.PeF.S,TF.*«M,TCIC) ► 

ViSI(n«TwnVARtNPkFS»MFMP,PLO(,,t ( H , P^^F S # 1 F , y l C J F 

(;u TD t7^* F 

cpicnscpsm F 

Tcu'ii n )stroNsa ) F 

visi (navissr n F 

17P DlFFCn>Sr>,HaUF-H*CTn)F*l,S)*^llN./(PF fM)*F l*FSl F 

IHP CONTINUE F 

00 IS|,nN f 

IF TCS(l)-l.*') F 

iNP continue F 

CP"'! iisrpsn ) F 

vis«(T>sviss(n F 

TcoN‘4(i>sTrnNscn ► 

un T j ?i*i F 

C^'^nTinuf f 

cp-^t nscpi f lucsf n*(cps( n*cpu i n f 

UU( 1)S( ( t .HtSURT ( Vissn )/VTSI F 

J H4*iU«T ( 1 ) T F 

1 ,/?! r n = ( n .vusuRT ( vTsi f 

I T r 1 /nmst ) f 

4S(nsfC5fnFrtMn/(-**^3Frs(i)»(«" f 

»» 1 1 nst 1 ) ► 

IF (X[( TT.FU,^\v1) (,.) T*) ► 

vrs^^n)aviss(i)/n,^^f, i^(n^*ni)/»s(in4vi'»irn/n,*^f-r?im*ni ► 
I S(l)/tl(l))) F 

TcUNM(il*TL0N5(n/(l,i4ft,viNS4(*l^r !)•(!! m/>«b(n)lFTCNNl(I)/(l, F 
1 P♦U‘4b5*c.^l (n*(xs(n/xKnn f 

continue f 

IF (IPE»F.F'^tn GO T*1 ^^U F 

call holier (M,Pt(H)#2#Tf n.^#tNT,WHOSn)»0) F 

RHOWf IlBWHUStt) F 

GO TO F 

HHUM(I) XPE (n)#««hs*hHI / (WUNI V* F 

continue f 

CONTlNllE f 

C F 

C •— •••COHPUTATtON OF FLO* PWnPFPTiFS AT Tr^F w, L , tOoF F 

C F 

ISN F 


eu 

ib 

^7 

ill 

SI 

i<e 

ii 

54 

ib 

ib 

57 

5h 

49 

'lU 

41 

4^ 

w 4 

44 

4b 

40 

u 7 

at 

u9 

s ^ 
bl 

b4 

b4 

bn 

-»/ 

s« 

bN 

Nl 

o 5 

r / 
n 

;/ 

7 \ 

7tf 

75 

7a 

7 b 

7n 

7^ 

7n 

N,, 

\ 

M.i 



o r» o 


TcoNMn>iTcoN8(n f 

VlSM(l)aVlSS(n f 

cpHn)*cps(i) f 

IF nP€HF,eg*n co tp f 

CALL ^OLiew (H,PE(^)#2*T(n,2»FNY«MHUS(I)»6) F 

MH11M(I)«RmOS(I) F 

UU TO 280 F 

25P PMOMI X )»Pe (M)ifci(M8/(BUNl v*32, ITaaTd ) ) f 

280 continue r 

BOJ-UE»PHOt (H)*VI8E(M) F 

OU 2 70 !al,fi F 

CWF(l)»NHOHfn*VlSMtn/»OHUF F 

;>CNO(n»viSM(ll/(WH(jH(i)#niFFCU)) f 

pwNO(i)»vis><n)*cpwn)/TcuNM(i} f 

270 CONTINUE F 

F 

COMPUTATION OF TmE HAT»IX EL^'^feNTS F 

F 

OU 2B0 Is2fNN F 

ALFTaALF#f 1 ,0*FLOAT(I»1)*0£LT) F 

Add# n*-rcMF(n*(ALFT*A2)/f2,0*finLT*A2n^(ALM/(a,%i*ntLn)*( F 

1 ALF*CPF(I)-FLC(n-ALFTA(CPFn^l )*CWF(1.1 ) W C 2 , 0*OEL T ) -2 , 5 ( H ) * F 

2 OFLCoSm)) F 

A(I»U2)s0,0 F 

Adrl#3)a0,fci F 

ACI,2f l)«CPF(I)*(UFfM)AA2)*fAlFT**2WfFri4n-F ( I • 1) ) / f CPH ( H * T T F 

I E(n:*«,0*(OFLT#*2)*777.O) F 

A( I »2#2)»»(t»Fn ) A ( alFT** 2) / r 2,ii#(0FL T A*2)#PNN0Cn 1 ♦ ( ALF T / ( 4 , 0 * F 

1 DELT))ACALFACPFd)/P«Nnf p«/AlFT/CPMrniA((CKFd4n<CPM(Itn/Pk F 

2 NO(l^n«CWF l*l)*CPM(I-1 )/■ *^N0( I •!)>/( 2,0 AtULT))«FLC(n-2.0*S(M F 

3 )*oFLCi)S(n*cpF(nA(cPS(i)»cPim)«ALFT*(Cb(i*n«ci>ci-n)/(CP^'( f 

8 n*SCN!J(I)A2,0Al)ELT)) > F 

A d #2# 31a( ALFT**2)*CPF ( T ) * f f PS ( P »C P I d ) ) • ( T ht I A ( 1 ♦ 1) -THfc 1 A ( I - 1 F 
I ) )/(CPM(n*SCNOd )*(B,oi*fOFLT** 2 ) ) F 

AfI,3#nsO,0 F 

Ad,J,2)»P,tl F 

All# J#3>»-(^ALFT**2)*C«F(n/(2,0*fnFLT•*^)*SC^0f n ) ♦ ( ALF T / ( a , 0 * F 

1 UtLT))*(ALF*CPF(I)/SCNOm*FLCd ) -2 , 0*S (« ) *'^FLCnS ( M • AlF T * ! CNF ( F 

2 Ul )/SCNOd*n*CWF (I-n/SCNTit T-1 ) W(2,0*iitLT))l F 

HCI , 1 , p»f alFT** 2)*CPF (n/(0fcLT#*2 WHt T Af M)*Fd )f2,.US(^) aF fl )/ 

1 OS(M) f, 

Hd#l,2)*vl,fr^ F 

B C I d , 3^ *0,0 F 

Hd #2» I ^*Hf TA (H)A( Ufc (M) a* 21*«H0F (M W(CPM( I ) A T1 E f «) ♦NHflNC P *77 7, F 

I P) F 

H( I #2,2)»t aLF Ta*2) aCPF (n/(PRN0(nAtDF.LTAA2nA2,0ASCM)AF(n/OS( F 

1 M) F 

Hd#2#3>*0,0 F 

Hd , 3, P«0,o F 

»d#5#2^*0,0 r 

ftd»3#3la(ALFTA*2)ArPF(n/(sCNUd)*frFLTA*2)lA2,0*S(MlAF (P/OS( F 

I M) F 

Cd # I r 1 )i*Ad # 1 f I )»(ALF T*A2)*rWF ( I ) / ( DfiL T A *2 ) F 

Cd»l#2>«O,0 F 

Cd#l#3)«O,0 F 

cd#2#n««Ad#2#n f 

Cd#2#2)«»A(p2#2)-(ALFTA*21*CPFCn/{PPN0(nA(i;tLTAA2) ) F 

Cd#2#3)*-Afl,2,3) F 

Cd#i»1 ^*0,0 F 

C(I#3,2)«0,m F 

Cd#3»3)*-A(I#3#3)-( ALFtAA2>*C«F(I )/f.-,CNJi(nA»OELT aa2) ) F 

Od#l,t)*CALFT*A2)*CPF(I)*(F(I + n«2,AF(I)iFd.p)/C2,UA(OF.LTAA2 F 


IM7 

IBB 

189 

190 

191 

192 

193 
t9U 

195 

196 

197 
19B 
\99 
2.V 
2«M 
202 

223 

224 

2 OS 
226 
2tf 7 
226 
229 
210 
211 
212 
213 
2P 
?tS 
216 

217 

218 
£:19 

220 

221 
222 
223 
22i 

22s 
226 
22 7 
22>^ 
229 
2 30 

231 

232 
253 

234 

235 

236 

237 

238 
23« 
2U0 

291 

292 

293 

p44 

295 

296 
2«7 


ORIGINAL 3‘AGB 1 ® 
OF POOR QUALITY 


o o r» 


1 )UtALFT*(F(Un*F(l-1 >)/(4i,v4*OfctT)1*(FLC(l)tALFT*(C»<^ (U\)*C«F 

2 (I.n )/(2,0*OELT)-AlF*C»F(n42,'^*S(w)*f>FtCnS(nUBETAC^)*WHOE(H 

) ) /KH0*«(n^2,k»*5CH)#(F ( n** 2 )/ns(M) 

l)(I#2f n«ULFT**2)*C»F(T)*(TMfcU(Ut )-2,i^*THfcT4n)iTHETAa-n)/ 

1 (2,fcJ*(DELT«*2)*PWM)(l)UCAl>T*(TNETi(nt)*TMfeTAn*n)/(<4,ir<*OElT 

2 ))*(2,K**S(M)iDFCClJS(n^FtC(n-ALF*CKFin/PW^'U(ntALFT*((CWFf ifl 

i )#cPH(wn/PkNo(itn-c«F(i-t )*cPHu-n/pwNu(i-n)/(2«v'40ELT)i/c 

<* PM(i))^^,tS*8(M)4F(n*THP.TA(n/0S(M) 

o(Ui,n»(AuFT#*2wc«F (l)♦(cs(I♦^•?•M•cs(I)♦cs^-m/(2,'/»•(0feu 

1 T**2)*8CN0(n)^(4lFT*(CS(Un*CS(l-1))/f«,-<*0FLTn*(FLC(IU2,k^* 

2 S(»il*nFLCnsri)*ALF*CPF (n/ 8 CNH(I)tALFT*(CWF (!♦! ) / SC^O ( 1 ♦ 1 ) *C«F ( 

3 l«n/SCMO(T •!))/( 2, »<»*lJELT>)42,W*S(^i)*F(T)*CSn)/i)S(*^) 

281^ CUnTI'^UE 

ENS«(^OVl^(M)*SCNO(n*5tJWt (2,*S(M))/(41 F PhOm ( n * V T Sh ( 1) ♦UM ^ W f ro 
IS(H)**hk) )i#(nFCT^.S*DELT*nfLT*n ,-(SCNU( n/C^^F ( U ) * ( C«F f 2 1 /SCNO f 2 
2)-C«F(n/8C'^O(n)/0ELTfk{WIN(M)*S»^»^T(2,^S(*^) ) *5CMM P / ( A^F *RHOh ( n 
3AVlSM(n*Ufc(H)*CwOSC<>**'^»v>)>^ 

tVALUATlON OF THE COEFFKIF'^TS 

Hf t , 1 , 

H( 1 r l ,2)»2,V^ 

G(ui 

t,(t,2,naTMFTA-i(^) 

G(l# 

Oil ia2,NN 

DH Jaws 

00 29^ 

AH( 1 # Jf K )aw, 

1)0 L*i,\ 

Art(I^J#K)a4M(I»JjK)4A(l § I^L l“l #Lf^^ 

29*1 CONTISUF 

oi) ’ 1^0 jat#3 

DO lat , i 

‘HAHf 

J«C1 Cor*iTI lu£ 

AU JHAHf I,l,naHAMf I,2,2>aR*M(l,5,5)*HAHn, J,21*MA^n,2,5) 

A 0 JBahO, 2 , t ^ 3 -BAM( T , 2, 1 1 *i 3 Ari( 1 , 3 , 5 ) 4 HAM( U 2 , 5 )*HAH (1 , J, I ) 
AOjrtAHfI,5,l)a8AH(I,2,n*HAH(l, i,2)*eAH(t,P,2)*hAM(J,i,n 
AOjHAHfl , I ,2)a-BAHn , I ,2)*MAHf 1 , 5f iWBAM( I , 1 ,3)* s 4H( i ,i,21 
A 0 JH AH Cl, 2, 2) 88 ah ( I , 1 » n*HAHfl , < , 5)*HAHf 1 , t > $) * 9 A H ( I , i , 1) 
AOJOAmC 1, 3,?)3.8AHn , I , naHAhU ,3,2 WHAhU , 1 ,2) AMAHM .3, n 
AOJHAH(I,\,3)88AH(I,l,2^<iHAHn,2,3>*bAMCl, W3 ^**^Ah( 1,2,2) 
AOJriAHC I ,2, 3) a«H AMU ,1,n*HAhn,2,3WHAM(I,l,3)*t5AH(I,2,n 
AOJHAHCl,3,3)aHAHn,l,l1*HAwn,2,2)*HAMf I,2#n*HAH(l,t,2) 
DeTHAHaMAHn,l,nA8AH(I,2f2)*MAM(I,.A,3UHAHn,l,2)*PAMn,2,5lAB 
\ AH(I ,3, IWHAHCI, I ,3)*hAHf 1 , i,?)*HAHf I ,2, t )-HAH( f , 1 ,3) aHAHC I .4,1 

2 )aH4hCi,2,2)*8AH( 1,2,3)*8AH( 1 , 4 , 2 1 ah ah ( T , 1 , 1) *H AH ( 1 , 1 ,2) aha h( 1 , 

3 2,l)A8AHn,3,3) 
i >0 5H Jai,3 

i)f! 310 hb\,3 

8AHlNVCl,J,K>aA|)JBAMn,.|,K)/DETHAH 

3lP CONTINUE 

DO 328 Jal,i 
DO 528 Kal,3 


F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 

F 


2*48 

2U9 

2b8 

2Sl 

2b2 

2bi 

2ba 

2Sb 

2Sb 

2b7 

2 bM 

2bP 

288 

2 bl 

282 

2ft4 

289 

28S 

?e 6 

287 

288 
28 R 

270 

271 
2/2 
2/3 

27 a 

27S 

2/8 

2/7 

2/8 

2 /p 

2 h»^ 

281 

282 

284 

2 «a 

2 Hb 

28e> 

28 / 

28 h 

289 

290 

291 

292 

293 

294 
29 b 
29 h 
297 
29n 
294 
300 

301 

302 

303 

304 

305 
308 
3«7 

308 

309 


00 no L«U3 

320 cnNTlNUC 

00 330 JB1#3 
AGUij^ns^t 
00 330 L«l»3 

Aua#j*n»AC(i#j#n+A(i,j,L)*Gn-i,i#n 
330 CONTINUE 

00 340 J«W3 

OA6(TiJtnaO(UJ,l)«AG(I.Jin 

3«0 CONTINUE 

DO 350 J»l,3 

1)0 350 L«lf3 

Gti,wM)«G(T# J# n^«*N!Nvni j,u*nA6n,L# n 
350 continue 

360 CONTINUE 

c 

C solution for the wi vector 

c 

«i(N,unsi,0 

*»fN,2, n«THET AE(M) 

^(N#s,n*i,0 

E (N)«W(N, \p\) 

THFTA(N)«*^(N,2,n 
CS(N)a-»tNp3,n 
DO 3*>0 ;I«Unn 

i«N-n 

00 370 J«\,3 

i)0 370 L*U3 

370 continue 

f )0 3M0 jai r 3 

I I * 0 0 1 ) 

3H0 CUNTINUF 

THETA(I)a«(Ir^f n 
C8(I)s«(i,3,n 
390 CONTINUE 

CONV(,s(F C 2 )*F n )-F 2 r 2 )^F^f n )/fF (P)*M n ) 

PRINT 760# CONVG 

IF t*<L00P,LT,2) UO Tt ao0 

IF ( ABS(CUNVG) ,LF#0#0005> r,o Tn 

av40 DO 410 I«lrN 

F2(i)«Ff n 

THETA2n)3THETA(I) 

CS 2 (I)sCSm 

410 CONTINUE 

DO 420 lalrN 

F(nt0,5*fF?(I)AFl(I>) 
TWETA(I)»,b*(THETA2tn^THETAl (1)) 
T(n»TTEfhlATHETA(I) 
csn)s0.5*(c$2n)4csi (in 
420 CONTINUE 
430 CONTINUE 
C 

C EVALUATION OF THE STREAM Fu nC T T (kw r E CL2 

C 

no 440 1«1 rNN 

FF(n«F (n/( ALF*(l,0*FLnATn.n*OELT)) 


F 310 
F 311 
F 312 

F 3 :; 

F 314 
F 315 
F 316 
F 317 
F 31H 
F 319 
F 320 
F 321 
F 322 
F 323 
F 324 
F 325 
F 326 
F 327 
F 326 
F 329 
F 3 30 
F 331 
F 332 
F 333 
F 334 
F 3 35 
F 336 
537 
F 33f* 
F ’339 
F 34^ 
F 341 
F 342 
F 3«A 
F 344 
F 345 
F 346 
F 347 
F 3 46 
F 349 
F 350 
F 351 
F 352 
F 353 
F 354 
F 3S5 
F 356 
F 35 ; 
F '.4 
F I 9 

F 3h»' 
F 361 
F 362 
F 363 
F 364 
F 365 
F 366 
F 367 
F 366 
F 369 
F 370 
F 371 





0 £LSTn)»UELST(I-2)*r*CT0W*f (HMnE(M)/KHnM(i-i).F{I-2))/(l,>».rN( F 

1 i*2))+a,»)*('»Hi)e(H)/RM0Mn-n»F ( i«n i/n ,»)- yn( i-n f 

2 M{n»F(i))/n,»i-YM<i)n*oEi.T/ti,«)**iF) f 

IF (n+n,CT,SN) so Tt) SaH F 

l>EL8T(I + l>«nEL8Tn-n+F*CT()B«f (HMOE f M ) /UMOm ( I- n-F ( 1-1 n / (I ,*«-r F 

1 lY{I-n)*«,0t(9H0EfM1/HH0M(I J-F(in/n,*«-YM(I))F(«HOEfM)/HHO»*(I* F 

2 n«F(un)/n.k'»vN(iYnn*DELT/(3,ia**LFi f 

5«0 CONTINUE F 

O5T*«(M)inOST»(M)*0EL8I (NN) F 

Th*10m{ F 

THMUi'(2)»F*CTOB«r fF m-F Cn*F( 1 n«CF C^T-F (?)*F(^T )/U ,»i-l)ELTn*OFL F 
lT/(2,H*iLF1 F 

DO 8SO IB3iNN,2 F 

THHUM(I)BTMMOM(I«2)*FACtf)H*( rF (I -? ) -F( I -2 ) »F ( I 1) / (1 . Y N ( t -? ) F 

1 )4<j,«)*(F(i-n-F(i-i)*Fti-in/n..i-YNfi-inF(F(i)-F(n*F(i))/n, f 

2 ki-VN(n n*OELT/(i.»'**LF) F 

.F ((Un.ST.NN) Go TO S52 F 

TM«uM(i*nnMMUM( i-n+F *cTOB*f(Fci«n-Ffi-n*F(i-in/(i.o-YNn« f 

1 n)»u.i*»fF{n-Fm*F tin/d.n-YN (m + (F fit n-F (1*1 )»F(un 1/(1, f 

2 0«YN(Un))*()ELT/H,»'*ALFl F 

sso continue f 

NSP*CtlY5t-N)/6 F 

IF (N8PfCE.lt,, 11 NSPACE»P F 

PBInT »v,i f 

('0 5o<) I«l ,NSP*CF F 

PWINT 7,10 F 

bhU CONTINUE F 

PHIM 77/, F 

PBInT 780 F 

PMINT 7B8, { I, YN(1 1 , Y ( 1 1 ,ET t( 11 ,F ( n , THfcTAd 1 ,CS(I ) , la1 ,N) F 

tPSd )tT"OVFH(NP9ES,NTFixP,PL0G,T d 1 ,PKEb, TFMP.CPSCl F 

YISS( n»Tt«nv*W(NPWES,NTEhP,PtOG,T ( 1),PF£8,TF0P|VSC1 F 

TC(1NS( ll»1-0Y*w(NPWEb,NTEMo,PLOG, T ( 1 1 , PNF. S » Tc mP , Tc SC ) F 

IF CAEY.EU.l ) UO TO S7o F 

CPU 1)sT«uv*R(npwES,nTF.mp,plOG,T ( n,PH£S,TFOP,CPIC) F 

TCONI ( n»T«0V4R(NP«ES,NfFNP,PLllG, T d ) , PBf S , T f oP , T C I C 1 F 

VISI (naTN0V4N(NPBES,NT£«lP,Pl.()r.,T ( 1 l,PPtS,TtoP,VlCl F 

(,U TO S80 F 

570 cpj ( n»cpsd ) F 

TCONI d )*TCONSd 1 F 

VISIdlaVISSd) F 

580 CONTINUE F 

CPMd JaCPI d UCSd )*(CPSd 1-CPI d ) 1 F 

IF (CSdl-l.O) 500,600,800 F 

500 XS'na(C3dl*i«Hl)/(wMS*CSdl*(«0I-«HSll F 

Xldlal.O-XSfll F 

IF (Xldl.Efi.O.Ol GO TO 600 F 

Gl2d )S( ( 1 .WfSUMT(YISSd 1/VISI ( 111»((NF I/t<*lSlP*0,2511»*21/(2,B280« F 
tSURId,*»MS/NHlll F 

G21 d )a( d ,n + suPTf YI8I d 1/vISSd 1 ) * ( ( wnS/nmI ) •*,). 25 1 ) **2 ) / ( 2. »2«o* F 
ISQHTd.FNMl/OMS)) F 

ViSMd laOiSSd )/d ,0«.&12d )**Jd)/ASdl)FYTSI( 1 1/d ,0 + G2*. d )»xSdl F 

i/*id)i F 

TCONMd )aTCONS( 1 l/d.O*l .065*G12( n««ld)/*Sd}l*TC 0 Nldl/d, 0 +l ,0 F 
i65*G2i(i)*Y8(n/*nni f 

GO TO 6H F 

600 V1SH( naVlSSCjl F 

TCnNMd)»TC(JN8d 1 F 

610 continue F 

T*UbbhoM( lUVISNd )»UE(M1*UE( 0 )»(Ni)S(m1**kK)*AlF* (F (21-F ( I ) J/(l)F.LT F 

1»8Q«T (2,»8(M1)) F 

CFaT*U/{0,5*0MnE(Ml*UE(Ml*uE(ol ) F 


BiB 

b35 

B36 

B17 

B38 

BiO 

B40 

BBl 

bb2 

BBi 

B4B 

UB5 

BBe 

Bb7 

auH 

aao 

aSt 

«S3 

abu 

ass 

aS7 

aStt 

aS9 

aferl 

I 

af}'4 

4r>S 

V 

4*^7 

UbO 

u 

u7l 
472? 
a/i 
u7u 
a7S 
a 7b 
a 77 
a 7h 
a7«4 

I 

46i 

a«u 

aHb 

a4b 

aS7 

use 

a«s 

a9 1 
U9i 
U93 
u9a 


Ht4T«»TC0NMn)*TTf H * ( N{)5 f H > * *»^^ ) * ALF* { Tm£ T A ( ? ) •TMt ^ 

JTAn) )/( 0 ELT* 8 QRT( 2 ,t 8 (M) )1 F 

QOOTs-HEAT F 

CPe«TrtOVA»(NP»ES^NT 6 MP,RLL".. Tt (M) ,P»ES, TEMFfCFSC) F 

rCONC«T*<OvAP(^iPRE 8 ^NT£MP,^ & r ^ i m ) , PKf 8 # T£HP , T C SC ) F 

PpNO(N)«VI 8 £(M>*CF£/TCUNt F 

«£C»SfJRT(PPNO(N )1 F 

TREC«M€C*(TTF(*^)»r£(M) WTE(*<) F 

5TN0sgtM)T/(PM0t(H)*U£(M)^CPM(N)*(TRFC*TMET*wf^J*TT£C^) ) ) F 

IF (AHS(ThETA»i#-U,LT, 1 ,E* 8 ) GU TU b^fc) F 

MC0EF«»00T/(TTECO*n.-THeT4iN(M) ) ) F 

r.O TO F 

hC0£F»99P^ F 

63V4 mRaTIO«HCUEF/HT«£F f 

PHfNT TAu»CF,MCOFF F 

PRINT HH, TREC#«0^^» 8TNf>,HP4no F 

PRINT DFLST(NN),THHf)H(NN),OSTAR(M) F 

DU bR0 I»l#NN F 

FF(I)aF(I)/(ALF*(l,fc1«FL04T (I*n*0£UT) ) F 

bRi^ continue F 

IF ( CO TO bSfci F 

FLC(l)*^^,5*(FLCw(N)4FLCi«(Mtn) F 

FLC(2)aFLCf 1)4(FM )/ALFtF (2 ) / ( 4LF* ( l , kJ*l)F L T ) )1 tQEL T/2 , ^ F 

00 65W Ia3fNN,2 F 

FLC(I)»FLCCl»2)^(FF(I»2>f«,0*FF (I-1UFF n))*0ELT/5,rf F 

IF f Cl 4 l),r,T,NN) GO TO bS^^ F 

FtC( l+n»FLCn-t J^^'FFfl*! ) + u,^t*FF HWFFn^l ) I*0ELT/3,^ F 

6 SP continue f 

no bOM 1*1,NN F 

OFLCoSCn»(FLCn)-FLri(I))/nbfH) F 

bb^ continue F 

FLC(N)aid.v 1 F 

i>0 b70 I*1,N f 

FI (i)«Fri) F 

TkctaI fnaTHETAn) ^ 

csi(n«cs(n f 

FLClv n-FLcU) F 

870 continue F 

F^miFUj) F 

F2f2)«Ft (2) F 

rtSN^I F 

GO Tl» \*y^ F 

CONTINUE F 

RF-TUrN F 

C F 

bPO FORR 4 T (Inn F 

7P^< format (/) F 

7p^ format (52X, 52HSTATIUN M*l, THE INITIAL PROFILE) F 

7^^ FURM4T (43*# RHF( I)#8X#MHTmET 4( n, 8 X,SHCSn ) ) F 

73 ^ format .9X, U5»3E1^.S) F 

74lrf format c//4J«X,flHTAu»rEl2.S, |^ix,3mCF*,E J2.51 F 

75P format (//3PX#5HT«EC«f n2,S,9X,5HUl)uT».El2,S,9X#bHSTNU»,F li?,')//) F 

7bM FORMM (1HX,6HC0NVC«,F J2,5) F 

770 format (50X,2HM«,I5#5X,5MAfM)«.FlU,b#/) F 

780 Format ( iSX, IhI ,2X,5mVN(1) ,2X, RMV ( I ) , 3x# 6HFTA ( n ,5 x,RHF ( P ,bX, BfiTh f 

lEtA(l)#SX^SHCSn)) F 

790 FORMAT (31X, IS, 2X,F«>. 3, 2F8.S,3F12,b) F 

800 FORMAT t//30X,RHTAU»,EI 2,S,bX, 3HCFt,Ei?,S,7X,bHHC0EF«,fc 12# S) F 

810 format (//30X,5HTREC»,tl?.5,SX,SHQnOTB,El2,S,5X,5HSTNU»,El2#S,5X#7 F 

lHHRATI0»,fc'I2.5//) F 

820 format (3nx#bH0FI.ST*,Fl2,S,8x,8HTHM0M«,F 12,S,bX,8Hf)ST4R»,E 12.S//) F 

8 S 0 format riRXP5,2I5X,FlS,8)) F 


U9t 

«97 

U98 

499 

500 

b01 

5^2 

521 
S<3« 
S0S 

S&*9 

SI 1 
S)^ 
SI ' 
SIR 
SIS 

SIM 

SI f 
SIM 
SIS 
S?0 

522 

523 
S 2 u 
S2S 
S2b 
S27 

S 2 H 

S29 

S31 

Si2 

S53 

S3a 

SiS 

S3m 

S37 

S3m 

S3^ 

SRk' 
Si* 1 
S**2 
S-»3 

%U4 

SRS 

SRb 

SRT 

s«» 

S49 

Sbw 

551 

552 
SSI 
SSu 

sss 

SSb 

SS7 


origin 

OF BOOR QOAUA* 






function TNOVAQdrVMX, JTyMMf VALX» VAtV#k* Y,FN) 

DIMENSION Y( JTVNX)^ XtlTVNx), FN(tTVNA«JTVMk) 

00 10 J«1»JTVMX 
IF (y(J)»VAtY) 

10 CONTINUE 

20 1>U SO Z«l#lTVNJt 

IF (X(l)»YAtX) io»atSSo 
30 CONTINUE 

FNMNSFNCIfJ) 

THUVARaf NHN 
RETURN 

50 IP1«I 
l«l»i 

TMOVARaFN(dJ)^(FNCIPt,J)*FN(l«.n)«(VALX»xrl))/(k(tPl)-xri)) 
RETURN 
60 J«J*1 

00 70 I«!#ITVmX 

IF UfP.VALX) 70,60fR0 
70 COnTInUH 
6i4 FNHNsFN(dJ) 

GO TO iMd 
0-^ IPl»I 
l»I-l 

fnhnsfkc 1 , j)4(FN(iPif j)*fn(i« jn*(vALx»k(in/(jf (ipn«%( in 

1 JUxJ^l 

00 110 isnnvHx 

IF (xm-vALo 13^^ 

1 10 CONTiNUt 
120 FNMX«F^{1 ,JiI) 

GO TO la'-* 

150 IP1«I 
I«I-1 

^NHXSFN(l,JU)^(FN(lPlfUl))«FN(nju))*(VAtX«X(n}/(XriPl)«X(I)) 
140 CONTINUE 

Tll•0VAReFNMN4(FN^:X•FNN^J)a(vAl,V-y(Jn/(Y(Jl'1•Y(J) ) 
return 


G 1 

G 2 

6 5 

G 4 

G S 

G 6 

U 7 

C « 

G 4 

G U> 

G 1 1 

G 12 

G 15 

G 14 

G IS 

G 16 

G 17 

C; 16 

G 10 

G 20 

G 21 

G 22 

G 25 

G 24 

G 2S 

U 20 

b 27 

G 20 

G 20 

G 50 

G 51 

b 52 

b 5 5 
G 54 

b IS 

b io 

b 5 7 
U i 0 


C 


END 



SUBMUUTIsI MOttEM (HpP«NOPT,T» 

C 

C ^lOPTsi) LUOK UP PPOPS BASED UN P AND H 

c nopt«i look up pnops based on p and S 

C NUPT«2 LOOK UP PPUPS RASED UN P AnO T 

C NOPTs) LOO«i UP PROPS BASED ON h AND S 

C 

dimension FtP(SS#2*^)» TT(53,2M), ZTr3S#2B)« GAMt(55,2B) 

1# ENTRO(i3,2D); FLPO(bS«*)« H20(bbB), TT0(6 nB), ZTUfbbB)* GAmE0(»«6B 
i)f ENTROlKbbP)# FLPZf?h')» HTHt(3S)» ENTUDV f S I) , FLPV(2«»»3n 

equivalence (FLPO^FlP), (NZOfM?), CTTu,TI), CZTO,2T). CGAMEO.GAME) 

J * (ENTwoOfENTRO) 

DATA PO,CPOR,MOfSO**#C,W,ALC.PTi1,CP/Zl 1 b , » 3, «« 1 5«, It 7 . 2 J, b r 5^. i 

l»Si,SS,2.3P2SRS»55,7VS,,21R6b/ 
iDEAlsbl 

2*1, 

G4MH4B|,a 

IF (I2.EQ,53) go TU 3B 
OD li* K 3 1,22 
LL*i3*(K«! ) 

DU IB L«lfl3 

FLPO(LL^t)*FLP2(A> 

IF CK.tO.I ) GO TO 12 

H20(U4L)»MZUa) 

U* continue 

12b33 

I2*2B 

DG 22 I81,IZ 

MTHU(I)»H2(I,n 

no 22 Jsi ,.!z 

entmovcj, nsENTNon#j2-j*n 
FLPV( J, n*FLP(l f JZ-JaH 
22 continue 

32 II- CN0PT,trj,3) GU lU 24«ri 
PL»AL0Gr2(P/?l 1b,)^12. 

IF (NOPT-O A2f|a2f2l2 
*2 IF Ch,LT,122,) go Tu 112 

S2 Continue 

Call DInT 2 (H.H2r2»21 » PL »FLP2f l2fJZ#T,TI ,SS»fcMQO, 2 , Ga^ma, game ) 

IF {N0PT,eQ,2) Go TO BV 

IF ( Am 5( (S»SS)#12, •*( 4-IPLACF ) ) ,LT, 1 , ) GO TU 92 
IF (S«SS) b2,92,72 

62 pll»pl 

PL»(PL^PLu)/2, 

GO TO S2 
72 PlJsPL 

PU«(PLtPU)/2, 

GO TO S2 
R 2 SsSS 
92 continue 

IF (2-,2222Ef 3P) 122,|32,liw 
122 RHn*P/(32,2*S3.3S*2*T) 

GO TO J32 

112 IF (M,LT,2,) go T(3 132 
T 3M/CP 

RHosp/(r,#»#T) 

122 S3(CPUR*ALUGt2(H/HO)*ALOGt2(P/Pn1 )*ALF4S0R 
CO TO 332 
I3«» IFP*1 
inEAL«i 
go to 332 


H 

1 

H 

2 

H 

3 

H 

tt 

H 

S 

H 

b 

M 

7 

H 

b 

H 

9 

H 

12 

M 

I 1 

H 

12 

n 

13 

H 


H 

lb 

H 

lb 

H 

1 / 

H 

IR 

M 

19 

H 

20 

H 


H 

22 

h 

23 

H 

2a 

M 

2b 

H 

2b 

H 

27 

h 

2H 


29 

H 

3o 

H 

U 

M 

32 


55 

H 

5a 

H 

3b 

H 

ib 

n 

37 

H 

3b 

H 

39 

H 


H 

a] 

H 

42 

h 

as 

H 

a a 

H 

ub 

M 

ab 

H 

U7 

H 

<*b 

H 

49 

H 

bo 

M 

SI 

H 

S2 

H 

bi 

H 

b« 

H 

bb 

n 

bb 

H 

bT 

H 

bb 

H 

b9 

H 


n 

bl 




!<*(^ IF >5 CO TO 

no I5i<» J«2,J/ 

IF (PL-FLP2(jn 
15^ CONTiNur 
l6bl IF (S*ENTHU(2, J) ) 

170 IF (S«ENTHf)(i,J«i)) 

leo CAtt OInT? fS,€NtWU#M,M2,PL,FtP7# W,JZ,T,TT,Z,ZT,l,GA«1M4,(;ArtF) 

IF 

14i4 H»»TO*cPOk*lO,**( fS»Sn»4AttAAtnGi^(P/pn) )/fCPO«*ALE)) 

T«h/CP 

PhosP/ rc«p*T) 

CO TO 
?P0 IFPsl 
I0EAL«1 
CO TO iS^I 

2to IF (T.UT.41P,) Gn TO 

CALL OlwT^ (T,TT,H,H2,PL»FtP7,I2,JZ,2,7T,S,ENTM0,l,CAHM4,CAM€) 

IF (H,,22p2F*jk^) Ulid, 25^1^230 
220 IF fT.LT.fr*,) CO TO 23n 
h«cp*t 

WM0*P/(C44*T ) 

uo TO 

230 IFP*l 
10EAL*\ 

Cl) TO 330 

240 IF (h.lT, 100,1 GO m 300 

call DINTI r5,FNT»0V,PL »FLPV,h,mTHL#sI2» T2) 

IF (PL*.2222E4 30) 2So, 310, 3io 
250 IF tPL,r,|,)?,<,0HP7) GO TO 32o 
IPLACt=AUnuiK(S) 

Iso 

2co Ui4l 

IF (FLP2(n*Pt) 260,270, 

270 PL'f=FLP2(l4i ) 

PlLsFLP/d-l) 

Go Til 2^#o 
280 PlosFLPZCn 

PLL»FLP2( !•! ) 

2P0 PS10,**(PL-lo,)*21l6, 

GO TO ti'A 

300 IF Th,LT,0.) 01) TO 310 

4LPx(S 0><-S) /ALF4CPf>^**LnGlofH/MO) 

PSP0#10,**4LP 

iBH/CP 

pHa«p/(G*k*n 
CO TO 33.1 
3!;’ IFP*l 
lOF aLs) 

CU TO 330 
520 I0F4U*2 
330 >ifcTUPN 

OATa (FLPZCJ), J= 1,20) /6. 301 «3# 6, 69897, 7,0,7,30103,7,69897,8,0,8,30 
I 10 3# 8,69897# 9,0,9, 3o 1 03 , 9 , 69897 # 1 0 , 0 , 10, 30 1 0 3 # 1 0 , 6989 7 , 11,0,1 1,30! 
203# 1 1,69897, 12,0,12,30103, 12.69897/ 
data (TtO(L)#L«l #2W4)/o,,ai9,,e3^,# 1616,,2539,,3032,# 3478,, 3898,,a 

1025.. 4175..432,),,496h,,U600,,tt8H«,,SU8,,S730,,6238,#6565,,683b,,7 

2005.. 7 1bS., 7270,, 7 380,, 7693,, 79 1 1 ,,81 00,, 8262,, 84 10,, 8568,, 672 3,, 8 

3900.. 91 40.. 9450.. 0.. 4 19.. 834.. 161 6.. 2359.. 3032. #35! 7., 3907,, 4130,, 
44300,#a46w,,a6O0.,4750,,49a0,,S260,#5760,,63O0,,670S,,6984,,7l95,, 

57350.. 7490.. 7596.. 7927.. 8 1 72, #8 574,, 8550,, 87 12,, 8878,, 9054, #9250, , 
5 9490, #9760, ,0,, a 1 9, #814, , l6 1 6, , 2519, # 3i*32, , 5546, # 3964, , 4224, #4 390, 

7#4563,#4708,#AH65,#5060,#5346,#S77S,#6550,,6795,#7o92,#730b,#7495, 


H 

H 

H 

H 

H 

H 

H 

H 

M 

H 

H 

M 

H 

H 

M 

M 

H 

M 

H 

H 

r* 

H 


H 

H 


N 

H 

►* 

N 

H 

H 

M 

h 

H 

H 

H 

H 

M 

►< 


H 

H 

H 

H 

H 

H 

H 

M 

M 

H 

H 

H 

H 

H 

h 


62 

65 

64 

65 

66 
67 
66 

69 

70 

u 

72 

73 

74 

75 
7o 
7 f 
78 
/V 
80 
81 
82 
83 

6U 

86 

87 

86 

89 

90 

91 

92 
9 3 
Qu 

95 

96 
9 / 

98 

99 

1 V •’ 

h 2 
iri 
I la 

I 

*'7 

i -8 

t.^9 

1 1 
1 1 
12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 
23 



u>* 

\ih 

U7 

131 
l3^ 
I3i 
1 3a 
l3b 
l3h 
13/ 
1 38 

I ai^ 
1<*1 
1^2 
1^3 
] uu 
laS 
1 

1«7 
1 Mo 
1 

ihl 
IS^ 
iSi 
iSa 
ISS 
1S^ 
1ST 
1 Sh 

1 *>t6 

IM 

l«>3 
1 ^a 
1*»S 

I 

/ 

168 
1 6*^ 
1 7^ 
1 71 
172 
1 73 
I 74 
I7b 
1 7ft 
1 77 
1 7ft 
1 79 

I8fc? 

181 

182 

183 

18 a 





b 119*1.1 47, 1.1 7 1.1,193,1. 2Mb, 1.21 3.1 

7287. 1.348. 1.4 11.1. 475.1.539.1, bMtt,l 
«1.8,l,nirtl,l.Mlb,l,237,l.MbM,l,k108,l 

9212,1.2l9,l,229,1.241.t,2SS,t.2bb,l 
*S9T, 1, 6b 1,1,72S, 1,788, 1.851,1,9^9,7 
•188, 1.132,1.15b, 1.181, 1,197,1.288,1 
•277,1.337,1.398,1.461.1,528,1.589,1 
•1.8,1,812,1.828,1,852,1,876,1.188,1 
•214,1,223,1.234.1.248,1.259,1.273,1 
*644, 1 ,7M7, • .778,1.831,1,891,7*1.8,1 
*128,1,144,1.169,1.188.1.281,1.211,1 
•328,1.387,1.449,1.51 1,1,574,1.638.1 

•818,1.822,1.842,1.863,1,888,1.112,1 

U*T* (7Tn(L),LB258,582)/l,2l5, 1 ,226 
11.439,1,588,1,563.1.625,1,687,1.748 
21,837,l,858,l,88a,l.lM4,l,128,l.l52 
31.232.1,244,1.258,1.312,1.372,1.431 
4 1,798, 1,855,7*1.8, 1 , 888 , 1 . 1 5 , 1,852 
51,164,1,184,1,197,1,287,1,216,1,228 
61.484,1.545,1.688,1.667,1 ,727,1.787 
71,. '46, 1,866, 1.888, 1 ,1 1 8, 1,1 33, 1,1 54 

81.232.1 .245,1.297,1.353,1.412,1 ,472 
91.826,7*1,8, 1 ,885,1.889,1.822, 1,848 
*1.164,1.188,1.193, 1.284,1.215,1,227 
*1,522,1.581,1,648,1,699,1.757,1.812 
*1,852,1,872,1,892, 1,115,1 .135,1,154 
*1,231,1.281,1.337, 1,394,1.452,1 .51.' 
*7*1 ,8, 1 ,8w3, 1 ,885, 1,815,1.827,1 ,844 
•1.159,1,174,1,188,1,199,1,218,1,221 
*1 .552,1 .689, 1 .666, 1 .721,1.777,7*1.8 
• 1.855,1 .872, 1 .892, 1, 1 12, 1 , '29, 1 ,148 
•1 .261 ,1 .314, 1, 368, 1 .424. 1 ,488, 1 ,537 

I'AT* (7Tii(li ,l* 58 3,6681/1 ,i'8l , 1 ,.V3 
I 1.882, 1 . 1198,1 ,1 1 7, 1.1 36, 1 .1 52, 1 .1 65 

21.355.1.418.1.466.1 .522. 1.577. 1 , «32 

31.814.1.824.1.837.1.852.1, /69,1,8H5 
4 1.1 76, 1.1 89,1.235, 1,285, 1,3 3**, I. 398 
51, 712, 8*1. 8, 1,81, 1,1. 8^3,1,818,1,8 19 
61,187,1,122,1,136, 1,158,1,163,1 ,1 75 
71,468, 1 .634, 1,586, I ,639, 1 ,692,9*1 ,8 
61.858, 1,864, 1 .879, 1 ,193, 1 .188, 1 . 122 
91 . 386, 1 . 357, 1 ,418, 1 ,46 8, 1 .512, 1 .56 4 
*1.819, 1 ,826, 1 ,839, 1.858,1 ,265,1,876 
*1.167,1.234,1.262, 1.352, 1 .382,1.431 

n*I* (G»*'£0 (l),Ls1 ,212)/2»1 .4, 1 ,398 

11.1 14, 1 . 188, 1 . 188. 1. I 14, 1 .1 31 , 1 , 169 
2 1,1 2 1,1,118, 1,1 8 3, 1,899, 1,894.1, '92 
31.111,1.122,1.138,2*1.4,1.398,1.344 

41. 113, 1.113, 1,118, 1.131, 1,1 62, 1.212 

51.1 16, 1, 189, 1, IMS, 1.896, 1 ,89b, 1 ,297 
61,124,1 ,136,2*1,4,1,398, 1 .344,1.517 
71.116,1,121,1.132,1,159,1,2.13,1,258 

61. 1 14, 1,1 89, 1,1 81, 1.899, 1 ,1 88, 1,1 81 
91,139,2*1,4,1,598,1.344,1,317,1,295 
*1.124,1,154,1,156,1.195,1.258,1,261 
*1.114,1,185,1.182,1.183,1.184,1.186 
*2*1.4,1.398,1,344,1.317,1,295,1,236 
*1.157,1.154,1.166,1,237,1.261 ,1,216 
*1,118,1,187,1,187,1.188,1.1 18,1.1 12 
*1.398,1,344,1,317,1,295,1,244,1,1 7b 
•1.154,1,181,1.226,1,258,1.227,1,182 


883 . 1 , 

221 . 1 . 

667 . 1 . 

115 . 1 . 

263 . 1 . 

1 . 6 . 1 . 
216 . 1 . 

652 . 1 . 

126 . 1 . 

332 . 1 . 

811 . 1 . 
228 , 1 . 

699 . 1 , 

136 . 1 . 

1 . 237 , 
1 . 888 , 

1 . 172 , 
1 , 492 , 
1 . 852 , 

1.248, 

1 , 844 , 

1 . 173 , 
1 , 531 , 
1,856, 

1 . 238 , 
7 * 1 . 8 , 
1 . 172 , 
1 . 569 , 

1 . 163 , 
1.278, 
1,882, 

1.164, 
1 . 594 , 
1 .888, 
1 . 179 , 
I ,*>67, 
1 . 183 , 
1.445, 
1 . 838 , 
1 . 223 , 
1,881, 
1 . 136 , 
l.MS, 
1 ,898, 
1 .468, 
1,344, 
1 . 226 , 
1 . 893 , 
1 . 317 , 
1 . 262 , 
1,898, 
1 , 295 , 
1 . 255 , 
1 . 183 , 
1 . 227 , 
1 . 286 , 
1 . 188 , 

1.164, 

1 . 174 , 
1 , 116 , 
1 , 144 , 
1 . 156 , 


016,t,84H,l,U67,t.M94,l, 

N 


232.1,244,1.256,1,272,1. 

H 

le? 

752,1,795,1.858,1.918,6* 

H 


148,1.165,1.188,1,284,1. 

H 


343, 1,4U6,I.47H, 1,534,1. 

H 


813, 1.832, 1.856.1,881,1. 

H 


226,1.257,1,249,1.264,1, 

H 

X^i 

715,1,776,1,848,1 ,899,7* 

H 


152,1.175,1,193,1,205,1. 

H 

\9U 

343,1.455.1.516,1.582,1. 

H 


824,1.847,1,878,1,896,1, 

H 


231,1,242,1.256,1,269,1, 

H 

197 

761,1,622,1,888,7*1,8,1, 

H 

l9b 

168,1.188,1,196,1,286/ 

H 

199 

1.249,1.263.1.316.1.378, 

h 


1 .667,7*1.8,1,889,1,817, 

M 


1.198,1.281,1.211,1,228, 

H 


1.554,1.616,1.677,1.738, 

H 


1.876,1.898,1,121,1,144, 

H 


1.252,1,386,1.364,1,424, 

H 


7*1.8,1,887,1,812,1,826, 

H 


1,166,1,288,1,218,1.228, 

H 

227 

1,592,1,653,1,711,1,778, 

H 


1 ,888, 1 , 101 , 1.124,1.144, 

H 


1 .269,1,345,1 .483, 1 .462, 

H 

212 

1.884,1.887,1,818,1.834, 

H 

2U 

1.186,1,197,1 ,286,1 ,219, 

H 

2lb' 

1.627,1 .685,1,742,1,796, 

H 

215 

1 ,881 , 1 ,182, 1 , 122, 1 , 148, 

H 

21** 

1.124,1 ,379, 1 ,437, 1 ,u9u. 

H 

21b 

1 .>’*.4, 1 ,812, 1 ,1*23, 1 .837, 

H 

21b 

1. 177,1. 169, 1.28M, 1,213, 

H 

217 

1 .658, 1.785,1,759,7*1.8/ 

M 

21b 

1,8 14, 1,8 38, 1,846, 1,864, 

H 

219 

1.191,1.283,1,251,1,582, 

H 

22^1 

1.739,6*1,8,1 , 8 . 12 , 1 ,885, 


221 

1 . i 14, 1 . 1 36. 1 , 158, 1 . 164, 

H 

222 

1 .58. , 1 ,553, 1 ,687, 1 ,668, 

M 

225 

1 .844,1 ,859, 1 ,875, 1 ,891 , 

H 

22 a 

1.272,1.322,1.375,1,426, 

** 

22S 

1 ,887,1,815, 1 ,824, 1 ,837, 

H 

22 b 

1 .148,1.161,1,286,1,256, 

M 

227 

1,665, 18*1, ,1,883,1,81a. 

H 

22B 

1 . 184 , 1 , 110 , 1 , 126 , 1 . 141 , 

M 

229 

1.531,1,561,1.638/ 

H 

25** 

1,517, 1,295, 1 ,286, 1,143, 

M 

251 

1,264,1.238,1,188,1.143, 

M 

252 

1 .894,1 , 896, 1 .899, 1 ,184, 

H 

255 

1 ,295, 1,214.1,149, 1,1 19, 

H 


1 .246,1,198,1 , 152, 1 , 126, 

M 

25b 

1 ,188, 1 , 182. 1 , 187, 1 . 1 14, 

H 

25h 

1,228,1,154,1.124,1,117, 

M 

257 

1 .196, 1 ,158, 1 , 1 33, I , 121 , 

H 

25H 

1.185,1.118,1,116,1,126. 

H 

259 

1,168,1,138,1,121,1,128, 

N 


1.165,1.148,1 ,126,1 ,118, 

H 

291 

1 • 1 12, 1 . 1 18, 1 . 128, 1 .148, 

H 

292 

1 , 1 36, 1 , 126, 1 , 125, 1 , 129, 

M 

295 

1,148,1,134,1.125,1.128, 

H 

29 a 

1.122,1.138,1,142,2*1.4, 

H 

29 b 

1 , 1 33, 1 , 1 38, 1 , 1 33, I ,148, 

M 

290 

1 .148, 1 , 1 31 , 1 , 1 2e, 1 , 114, 

H 

297 




I 

r 

I 

I 

I 

I 

I 

I 

I 

I 

i 

i 

I 

I 

I 
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*i,m,i,tn,t 

. 1 1 ?. t 

. lli.l.U»>.I, 

12/, 1.125,!. 133,1 

, 1 

.S’?/, 

6 

246 

*1.344,1,317,1 

1 

,^S*^.l,l«5,1, 

151,1,158,1,134,1 

, t 57, 1 , lain 

.Ibb/ 

H 

249 

OaTa (GAHEua 

>,L=^1 

J.?7M/i , 1 7Q, 

1 .218, 1 .2b2, 1.236 

<1.191,1.165 

,1.14 

M 

2bi8 

17,1,136,1,131 

. t.l 1« 

.t,ll“.l,ll«. 

1.115,1.117,1,122 

<1.1^5,1,126 

,1.13 

H 

?5l 

26, 1,147,i>*1,4 

. t. 

.1.54U. 1,517, 

1.29b, 1 , 2 s®, 1,193 

,l.l01,t.l“6 

,1.14 

H 

2b2 

ll’, 1,14?, 1,146 


. 1.177, t.,»,1<J, 

I ,244, t ,244, 1,2/7 

< 1 . 1 77, 1 . 159 

,1.14 

H 

2b3 

46, 1,139, 1,1?4 


. 1 , 1 1**. 1 , 1 

1.122#1.124,1,128 

,1.155,1.1“! 

,1.15 

H 

254 

5k*, 2*1, 4, 1 , 39.1 


,1.5l7,t.<»“S, 

1 .263, 1 ,2/1 ,1.169 

,1.tS<>,1.l“6 

• 1,14 

H 

255 

67/ 






M 

256 

!>ATA fUAMflUL 

) ,Ls^^ 

7,a^a)/i,is^, 

1.16/,1 ,177,1,2/4 

,1.^57, 1.^“S 

,l.?l 

H 

257 

19,1, 189, l^lb9 

f t . 1 ‘’n 

, 1 ,1“S, 1 . 1.7M, 

1.124,1,123,1.124 

,1.120,1.126 

,1.13 

H 

256 

22,1,137,1,143 

.1,1*’“ 

,'>*1.“.l.5'».>. 

1.544,1. 517,1, 29b 

,1 ,2o7, l.?ln 

,i.n 

H 

259 

37,l,lb9,l,iS2 

.l.is^ 

,1,1 So, 1 ,101, 

1.176,1 ,2/8, 1 ,23/ 

, 1 .2“2, 1.220 

, 1 .^k" 

H 

26ii‘ 

41,1, 179, l^lb4 

. 1 ,1^«? 

, 1 . 1 5-i. 1 . 

1. 127, 1,12«, 1,13* 

,1,1 52,1,150 

,1.1* 

M 

261 

S‘i,l,147,t,lSb 


. 1 ,5<J'', 1 , jaa. 

1.31 7, t ,29S, 1,271 

,1.225,1.166 

,1.16 

M 

262 

b6,l,lb.i,l,lS9 

.1,1'’! 

, 1 . lo'7, 1 , 10 1 . 

1 .197,1 ,221 , 1 ,23b 

, 1 .251 , 1 .21“ 

,1.19 

H 

263 

72,1, 174, |,ib2 

.1,1“' 

fl,l5S. 1,155, 

1 . 1 33, 1 , 1 3b, 1 , 1 36 

, 1 . 1 a<», 1 . l«o 

#1.15 

H 

264 

82,1,161,2*1,4 

. 1 . i**-' 

.1.5X4,1.517, 

1.29b, 1 .273,1.232 

,1.190,1.175 

• 1.1b 

h 

2o5 

97,1,1 64, 1 , Ibb 

.1.1^5 

, 1 . 10.’. 1 . 1 “1, 

t .21 7, 1 ,231 , 1.232 

, 1,222,1.2102 

,1.18 

H 

266 

*3,1.17^,1,147 

. l.iu-’ 

. 1 . 1 5», 1 , 1 50, 

1 . 1 4/, 1,142,1.146 

,1.lS^,l,l5o 

#1.16 

H 

2b7 

* 4 , 2 * 1 , 4 , 1 , 39 / 

.1,5““ 

.I.517,i.?9S, 

1 .274, 1 ,23®, 1,2/6 

,1,165,1.1 7^ 

#1.17 

H 

266 

*•^#1.171,1.174 

. 1.1««> 

.1.107,l.^lu, 

1.228, 1,232, 1,22b 

,l,2I0,l.l<»^ 

#1.17 

H 

269 

*9^1,lS3,l.iaS 

.l.t“i 

, 1 . 1U5, 1 . las/ 




y* 

2 7/ 

DATA (GAH^nfL 

1 , t=42S, b3b ) / 1 , 1 47, 

1 . 1 b, 1 . 1 54, 1 , 1 bk‘. 

1,1 06, 2* 1,4, 

1 .39fc' 

b 

271 

l#l. 344 , 1,317, 

1.29b, 

1 ,.>7S, 1 ,^aa, 1 

,218, 1 , 1 9S, 1 . 142, 

1.176,1.176, 

1.183 

b 

272 

2#l,l9n, 1,2/1^, 

1,21 3, 

\,iiS. 1 ,?5I , 1 

•228,i.21M,i,2t4, 

1.191,1.102, 

1 .153 

b 

273 

3, l,lS^,l.lSii, 

1.1b/# 

l,lS5,».lSo.l 

. 1 bvi, 1 . 1 h/. 1 . 1 7 3. 

2*1 .“, 1 . 59kt. 

1 . 344 


274 

4#1, 317, 1,29b, 

1 .27b, 

1 ..’'‘7, 1 .^.> 0, 1 

,2‘'4, 1,189.1,184, 

1 . 164, 1 . 10«, 

1,194 

H 

27b 

*># 1 .2*^2, 1 ,21 \p 

1,224, 

l,^5tJ, 1,^5 *1,1 

.224,1.212, 1 .22/# 

1.109,1.109, 

l.lSb 

H 

27b 

6,l,l5b,l.lSb, 

1.156, 

1 . 101 , 1 ,1oS, 1 

, 1 71 , 1 . 1 77,2* 1 ,4, 

1,59/, 1.5““, 

1,317 

n 

277 

7r 1 .29b, 1 .27S, 

1,258# 

1 .^5^. 1 ,ai 1 . 1 

.196,1.19,1.1.191/5, 

1.19i,1 ,1«6, 

1 

H 

2/b 

8,1.214,1,224, 

1 . 23 /# 

' ,?5l , 1 .?^7. 1 

.219,1.2/8,1.176, 

1 . lOS, 1 . toi , 

1 , lb6^ 

b 

279 

9,l.lbl,l.lb3# 

t , 1 b6 , 

1 , 109. 1 . 1 7S, 1 

.182.2*1.4,1,39/, 

1 .5““, 1 ,517, 

1 ,/9b 

H 

2«/ 

1 ,27b, 1 ,2«i3# 

1.25^. 

1 ,?l<». 1 'S, 1 

.199,^.104.1.199, 

t , 2 ^ 4 , 1 , <> 1 6 , 

1 ,/1 #> 

H 

281 

.22^, 1.25/# 

1 ,25<?. 

1 ,.>51 , 1 ,/^o, 1 

,2lb, ) . 1 8b, 1 .174, 

l.|o9,l.l'’7. 

1 . Ib8 

M 

282 

*f 1 , 1 7/, 1,173# 

1,177, 

1 . 1 " 1 . 1 . 1 0 7 , ^ 

* 1 .9, 1 . 59 1 , 34a, 

1.51 7,1 .^45. 

1 .^7b 

7“ 

2b3 

*• 1 .2b4, 1 .24/# 

1 

1 ,?1 ?. 1 .,?.-s. 1 

.2'^3, 1 .224, 1 ,2. ‘8, 

t .21 1 .22/, 

1 .^?5 

H 

?H4 

*#l,23t #1,235# 

1,^5^. 

1 ,^/0 , 1 ,.>.11 , 1 

. 194, 1 . 1 8 1 , 1 . 1 7b, 

1 . 1 75, 1 . 1 7«, 

1 , 1 7b 

H 

/8b 

*rl,W9,i,iA2, 

1 

1.1 .^»1 1 

. 39 . 1 , 1 ,3a4. 1.317, 

1 .290, 1.270, 

1 .#?bb 

b 

2**o 

*# 1 ,241 , 1,22b# 

t.^l7. 

1 .?! 1 , I ..>■■ 0 , 1 

.2vw, 1.212, 1 .218, 

I .22 5, 1 .220, 

1.229 

79 

28 7 

*#1.233,1,252, 


1 .??■*, 1 I 

.I8b,t,l*“2,l.l 79, 

1 , 1 «•’, 1 . 162, 

1,184 

b 

28b 

*#1.167,1 ,191, 

1 , 1'lo, 

<’•1. a, 1 . 59 * 1,1 

. 544, 1 .31 7. 1 .29b, 

1 .270. 1 .2SS, 

1,242 

b 

P89 

*/ 






H 

29fc^ 

oaTa (bA^ff'ri. 

) ,l Zb37,#)b5 )/l 

1 .223, 1 .22, 1 .214, 

1 . 2 1 6 . 1 . 2 1 9 , 

l./?3 

H 

29J 

l#l, 228, 1.228, 

1 ,^5M. 

1 .235, I ,2 3 5, I 

.229, 1 .227. 1.214, 

1.196,1,191, 

1,187 

H 

292 

2,1.186,1.1^9, 

1,1<<?. 

1.198,I^IQ«,1 

.2/2/ 



M 

d93 

OATA (eNT9hn(Ll#L*l 

,2»-i4)/23,bM,3 

1 fSb, 3^,97, 3b, 43 , 

57.42, 44,. 1 , 

39,86 

b 

294 

1 ,4/, 87, 41 ,78, 

42. bS, 

43,b''*,4a, 33,4b, lb, 45 , 94 , 4b, b/, 

“7, 55, “7, 94 , 

46,5/ 

H 

29b 

2#49,H4,49,SM, 

S/,18, 

S.’.S9,St 

3,^'‘b,S4,8«.bb, Tift, 

S“ , 40, o2, to. 

bl ,93 

M 

2?n 

3#b3.bl ,bb,22. 

bb, bS , 

00.9,1, 7o, 5 

/.64 , 3 5, »fb. iS. b 1 , 

57,//, 56, ,9, 

38, 9o 

b 
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4 A,.^l,ah^ 4 j,i 48 . 9 ^,a 9 ,a 3 ,a 9 , 93 ,'i 

W 4 -\ 99 ,ai , 79 ,a 4 >,S 7 ,a 

a 7 , lo,a 7 , 5 bf 


.•^H,au,rta,aS,Sb» 4 ^,i?^,ab,ei,a 7,30 
.^a,S 3 .Sft,Sb, 7 ^,bv^, 3 b 

.9*5, 3^, 3ft. ia,b^,3h, il , 37.4i^,3e,25 

, 3 ^,uti,Pl, 444 . 75 » 9 S,av^, 4 S, 99 ,ab,S 3 

,sa,3^#ss,99,sr,64,s<»,?<* 


t)*1,H^,6^,33^h3,ft«#bb, , 37,^9.Po, 3t ,f5e>, 3u, U,3S,bU ib,7il, 37.S6 

3 ^* 5 H, 39 ^aw’, 4 v\^ 4 ,iii ,.11 ,ui , 79 ,a?,S^,a 3 ,^^,a 3 , 9 ^>,uu,Sb, u'., ib,a^, 7 / 
Ub,2k1,4b,7vi,a7,l 7,u7,rtS,aH. 1 1^,4 9.94 ,S1 ,bS,S3. 51 ,SU,9b,b^,S?,b8,‘^6 
S<?,bU6l , n^6?,«i7,e>3,oq,^^^^4b,PM, 5U^ 3V , 75,3 5,?U iu, 7vi, 3S, 79, 3b,b3 
3^, S7, 3»^,a3, ,4 ',7 7,4 1 14, 4?, 6,^,4 3,44,4 4, ^2, 44,S7 

45,^9,4S,5S#ab,t^5,ab,4 9,4^,9 3,4M,7 3 ,Si% 3S,S1 ,9ft,5 5,ba,5S.^b.S6,Sb 
58,fcia,b<7,4b,bi^,Mb,h,>,,»rt, 19, 7 7,^7,i>5, 5i',v*b, 3^.S^,i4,^r'l , 3S,lid,3S,9b 


<? 9 rt 
<?99 
3 tUJ 
i^l 
iM 
3 ^ 3 
3/4 
3 .^b 
3 »^b 
3.1 r 

3 *"b 

31^9 


m. 



c 


*» j6,89.17,7i» H JJii 

• l8»«5,6j#«fc,*»4,«7,77,44,}8,S^.««,52,«9,SJ,97,55.«<' n JU 

• ,So.87»S8.?}«S9.«»^><»^.<>8#l9.**7,eb.96.^4,56,)1.8^,3}.3t>i4.a0/ M Jl? 

l)*t* (ENTfcOn(L).L«2l*5,<U*e)/|5,i7,36,l7, J7,i»t»17,«(#, J8.St,i9,i2, J9, m 313 

l9rf,aii,5 7,al,l7,«l,8i«««?,37,a?,<»t»,4 3.39,«3.»6,tta,32,««,7b,«S, J7,«b. h 3i« 
^88,».3,5^,^»9,9J,5l,‘•^»»5^.«7,5a,3^.5S,♦>6,•)7.8lr1,58.3b,59,68, H 3J5 

3^4»28.«S»3i«.9t»32.4i«,S3,99,34.3S»3S,^S>38.88#3«i.8$,3 7.S6,38.^‘i,38. h 31« 
«88,39,55»4**,l2,atf',ba,«j ,25,«t ,8H,42,i»'»tt2,7a,«i,lb,a3,61»9«,/M,uS, H 317 
5bl,«7,U,a8,82,S»i,^.»,51.<»i>,52.‘‘3*S«, 1 5,*»5,an,Sfc,68,37,9^, 1 7,«6,2S, m 3i« 
b3«»27,75»3‘’,21.31.7**,32,79,33,6«.3«.SS,3*>.38.1b.«2,36,8l,37,‘)H,38, h 319 
713,38,75»i9,33»39.9«»,at4,«3,4.4,95,«1,«^,«l,89,«2,i2,42,72,«l,13,*4«. h 32»i 
8 7.»,«6, 1 9,47,02, 49. H»'<,S»*,33.S1 ,ba, 52,97,S«.?<».5S,42,S6, bill, Ih, 77,24, m 321 
9b5,27,-lb,29,S2,3l ,•<» , 52. 1 8, 32,9b, 33, 88, 34 , 66, 55, 40, 3o,Vt8 , lb, 72, 37, h 322 
•37,37,97,38,54,39. t I .39.65,40,12,40,60,41 ,i«2,4l ,48,91 ,88,42,27,43, « 323 

»• J, 45, 25, 4b. 64, 47. 96, 49, 28, S0. 52, 51, 60, S3, 20, 54, 20, 55, 32, 15, MS, 25. m 324 
.7t, 26,14, 28. 60,30. ^9,31, 18, 32, .’3, 32, 97,53. 75, 34, 45,55,1 1 ,35,76, 3b, h 325 

• 3'-, 3b, 93, 57, 48. 38. 84, 38,54,59,02, 39,4 8,39,94,40,35,40, /5, 41 , 13,42, h 32b 

•60, 4 3, 49,45,5 3, 4b, 65,47, 8 7, 49, 06, 5-', 26, SI ,4S,52,SM,S3,b7,lS,lb,23, m J27 
•04 IS, 4'S, 27, 9i, 29, 4.1, 50, 49, 31, 34, 32. 26, 33. 04, 33. 74, 54. 41, iS, '.'ll, 35. m 326 
•bl,3t, 16, 56. 69, 37, 23,37,73, 38. 19,38, 66, 59. 10, 39. 58, 39, 89, 40,28, 41, n 529 
•70,43,04,44. 35,45,60,46,80,47,98,49, 10,50.27,51,35,52,40, 14,47,22, m 330 
•35,24,76,27,22,28,71,29,83,30,67,31,57,32.55,33,02,53,68, 34,28/ m 551 

n*T* (E8T«on(L),L»“09,6i2)/i4,85,35,40,3S,92,36,45,36,92,37.40,37, h 352 
184,38,27,38,68,39,05,39.40,40,81 ,42,15,43.36,44,56,45, 75,46,88,47, M 355 
290,49,07,50.15,51 .15, 13.55,21,45,23,84,26.30,27,79,28,88,29.75, 30, h 334 
364,31,42,32,10,32.72,33.30,33,87,34,39,34,90,35,40, 3S.«6, 3p, 32, 36. m 335 
4 75,37,17,37.57,37,95,38,30, 39,65. 40,93, 42. 1 0 , 4 3.2 7, 44 , 3a, 45. iS, 46, h 33b 
551 ,4 7,53, 46,53,49,50, 12,86,20. 74,23,15,25,61 ,27, 10 , 26, 19, 29,. 11 ,29, h ,3 7 

697.30.71.31.39.32.00. 32.60.33.15.33.68.34.14.54.62.35.10.35.54.35. h 338 
795,36. 35,36,77,37.1 1,37,44,36,77,40.02,41,16,42,27,43.32,44, 37,45. h 339 

642.46.38.47.36.48.28. 12.16.20.05.22.45.24.91 .26.46.. 27. 49. 28. 36.29, h 340 

925,3 1,03,30,68,31,30,31,86,32,4 1,32.9 3,33.40,3 3,86,54. 31, 3u, 73.35, h 34 1 
•14,35,53,35,91,36,27,36,62,37,90, 39,10,40,22,41 .31 ,42, 34 , 45 , 31 ,44, m 3“2 
•28,45,27,46.18,47.09, 11 , 25 , 1 9 , 1 3 , 2 1 , S4 , 24 , 00 , 25 , 49 , 26 , 58 , 2 7 , 44 , 2M , h 343 
•34,29,11,29,77,30. 37 , 30 , 95 , 31 ,46,31 ,96,32.41 . 32,86,53.30, 33, 72,34. 8 3«4 

•1 1 . 34 . 50, 34 . 84, 35, 19, 35,53, 36. 7 7, 37. 90, 30,99.40,00, 41.00,41 .95,42, h 345 
•87,43,79,44.63,45,49, 10.55, 18,44,20,84,23.50.24, 79,25,68,26, 75,27, 6 346 

•64,28.42,29.06,29.69,30,22, 30. 76, 3 1 .22, 31 ,68 . 32. 1 2 , 32, 54 , 32. 95, 3 3, m 3«7 

•31,33.68, 34,05, 34, 58, 34. 71 , 55,93, 37,,'8, 38,05, 39,08,40,0.1,4 - ,92,41 , o 340 
•8»r ,42.64,4 3.51,4 4. 30,9,86, 1 7, 76,20. I S, 22.61 ,24, 10,25. 1 9, 26,0«, 26,9 « 3“9 

•5,27.72,28,37,28,98,29.53, 30 , 04 , 30 , 50 , 30.95,31 .40, 31 . 80 , 52.20/ ►- 350 

OAT* (fNT»OII(l ),L«613, 6601/52, 55, 32, 91, 35 . 24,33,59, 33,90, 3S.. -9, 36, H iSt 
1 16, 37, 18, 36. 12, 39,. '4, 39 . 90, 40. 76, 4 1 . 5 7, 42, 37, 4 5, 18, 8. 95 , 16,8 3. 19.2 « 352 

24,21,70,2 3.19,24,28,25. 15, 26. '13, 26, 7 9, 27, 44, 28 , 0 3, 26. 59 , 29,1. -.29 , 5 M 35 3 

36. 30.00. 30.41. 30 . "0,31 .18,31,54,51 .89,32.23, 32.54, 32,84, 53,94, 35,0 h 35u 

41,35.96,36.90,57.77,56.57,59.33,40.12,40,92,41,68/ " 355 

P*T* (H20(L 1 ,l.»l , 53)/0. , 100. ,200, ,400, ,600. ,8.1.>, , 100.', , 1250, , 1500, " 356 

1. 1750.. 2 000.. 2250.. 2500.. 2750.. 3000 .. 3250.. 3500.. 3 75.1.. 4000.. 4250. " 357 

2. 4500., 4 750 ,, 500f'. , 600 . 1 . , 7000, , 6000 , ,9000, , 1 .3000. , 1 I 000 , , 1 2000 , , 1 3 -1 358 

3000 ., 1400!!,, 1S000,/ H 359 

M 

END 50 I 


nW POOB- 





n I 

c I 

C »)OUbLt iMf KPULATl^ii^ SDHWUOTlNt, U»l# TMfc SECDNO 1 sOE »^T I 

C variable is not constant "TTh the fHST, (THE RANGE Of The SFCOno I 

C IS hot TH| S^HE at each value Of THi ETPST.) 1 

C I 

OINtNSTON ZT(NL)» AT(Mt,NLT, YtlCHL^NtT, YTZ(ML»Nt), VTSCNt^NU I 

niNENSlON VT4(NL,NLT J 

IF (ZZ^ZTCn) 4v4,Ji^,lO I 

10 IF (IL.EO.n GO TO 2/ I 

IF (XX,uT,YTM,n) GU TU 4,fl I 

20 rn io i«i,NL t 

L»l I 

LL»l«l 1 

IF (2Z-?T(n> 1 

30 continue T 

40 vis,22?2Ea30 I 

RETURN 1 

SO RaTI(»»(ZZ- 2T(LL) T/<ZTa)-ZI(LL>T I 

iF(tL)lviO,l0O#l4O 1 

bO AO 70 J31,HL T 

L^*0 I 

LL^s 1-1 1 

IF (XY*»T(J,Ln BO, 40, 70 T 

70 continue T 

GO TO 40 I 

bv1 4ATI0£Un-XTaL^^*Ln/(xT(L^'»LT-xT(iLB,l )T I 

VlsYTULLN#LW«^Tin*{ VT I CL ^<#0 -V T I (LLB,l ) ) 1 

'^i?sVT2(LLN,| TfRATTn*C vT2CL*M.>-VT2(LLB,Ln 1 

yj3YT3f t LN#L)aRAT ( YT M L^’»L > • Y T 3 ( uLM , L I T I 

ytsYTi4(LLN,LT^MATlO*( VT4(L*^>L)*YT4(LL^#U T T 

return 1 

Yl«YTl 

v«>3YT2{L"’»Ll 1 

Y5»YT3(lb,LT ^ 

YtayTa(LM,Ll 1 

RETURN T 

UO *)0 no vUw*<L 1 

LHs.I I 

LL^aj-l 1 

IF rxx-xT(j,Ln i 

no continue 1 

GU to ao I 

120 yixVT I (L'I^LD-RAT lP*(vTnLN,LU-yT t (LM,Ln T 

y?*YT2(L4,LL)-WATlRA( yT2(LN,LLT-YT^(L«.L) ) I 

Y53YT3(LN,l LT-WAT IP*f vT 5(L^fLLT«YT3(LN,Ln 1 

tEsyT4(t4,Ll)»»ATlP*(VT4fLN,LL)-yT4(L«#U) I 
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Y9»vTJ(LL0»LW«*TIO«(yT3(lm,L)-yT5((.CM»L)) I 7n 

100 00 190 J»1.ML 1 HU 
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Y2sV6-9ATlP*(y6-y71 I <>6 
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ME Turn 

Pt' Vi«VM(LM^L) 

RF TURN 

UU 11*5 J = 1 ,HL 
L*^*J 

LLNsJ«i 

IE (XX«XT(J,Ln 
1 CONTINUE 
uo To an 
ia« TSSYT 1 fL'X,L> 
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